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ABSTRACT 


In  this  report,  we  describe  the  total  effort  at  The  Ohio  State  University  on  the  Project 
Analysis  and  Control  of  Interconnected  Structures,  initiated  on  September  1, 
1992,  sponsored  by  the  Air  Force  Office  of  Scientific  Research  (AFSC),  under  Contract 
F49620-92-J-0460. 

The  main  thrust  of  our  effort  has  been  to  consider  the  control  of  interconnected  rigid 
and  flexible  structures.  In  order  to  accomplish  this,  we  have  utilized  different  analysis 
approaches  including  partial  differential  equations  representing  simple  substructures 
and  finite  element  models,  radial  basis  function  neural  networks  and  wavelet  neural 
networks  representing  complex  structures.  Finally,  we  have  introduced  the  use  of  a 
multiresolutional  technique  with  wavelets.  As  far  as  control  techniques  are  concerned, 
we  have  concentrated  on  sliding  mode  control  with  a  detailed  analysis  of  sampling 
effects.  The  multiresolutional  wavelet  approach  which  we  concentrated  on  during  the 
last  year  has  naturally  led  to  a  multiresolutional  controller  design  in  the  domain. 
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IV 


OVERVIEW 


In  this  report,  we  describe  the  total  effort  at  The  Ohio  State  University  on  the 
Project  Analysis  and  Control  of  Interconnected  Structures  which  was  initiated 
on  September  1,  1992  and  sponsored  by  the  Air  Force  Office  of  Scientific  Research 
(AFSC)  under  Contract  F49620-92-J-0460. 

The  first  part  of  this  report  focuses  on  sliding  mode  control  for  sampled-data  sys¬ 
tems.  In  the  second  part  of  the  report,  we  concentrate  on  a  generalized  sliding  mode 
approach  to  vibration  suppression  in  flexible  structures.  Methods  for  block  control 
of  differential-difference  systems  and  a  sliding  mode  approach  to  control  of  dispersive 
flexible  structures  are  discussed.  Examples  of  each  method  are  also  presented.  The 
third  part  of  the  report  emphasizes  a  general  design  procedure  for  radial  basis  function 
and  wavelet  neural  networks.  Simulation  results  for  these  two  types  of  networks  based 
on  predicting  actuator  nonlinearities  occurring  in  a  flexible  structure.  The  last  part 
of  the  report  explores  uses  for  time- frequency  wavelets.  These  include  the  modeling 
and  control  of  structural  systems  and  subsystems.  Versatile  structure  model  based 
on  time-frequency  wavelets  are  developed.  A  multiresolutional  controller  which  min¬ 
imizes  the  sensitivity  function  for  a  SISO  flexible  structure  is  derived.  The  controller 
utilizes  the  structural  models  based  on  time-frequency  wavelets.  Implementation  is¬ 
sues  for  the  multiresolutional  controller  are  also  discussed. 

Contributors  to  this  three-year  research  effort  were: 

Students: 

Wu-Chung  Su 
Layne  Lenning 
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Shyamala  Raghunathan 

Post-doctoral  researchers: 
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A  number  of  publications  have  appeared  based  on  this  project.  They  are  listed  in  the 
Appendix. 

This  project  is  also  linked  to  AASERT  Project  F49620-92-J-0299  entitled  Integrated 
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Circuits  for  Distributed  Control-(AASERT  FY91).  The  major  portion  of  the  research  ! 

performed  by  the  students  supported  under  this  AASERT  Project  is  documented 
separately.  Some  of  the  developments  which  were  produced  in  the  AASERT  project  ! 

are  integral  to  the  research  reported  here  and,  therefore,  are  included  to  produce  a 
stand-alone  document. 


VI 


Contents 


1  INTRODUCTION  1 

2  IMPLEMENTATION  OF  VARIABLE  STRUCTURE  CONTROL 

FOR  SAMPLED-DATA  SYSTEMS  3 

2.1  Background  on  Sampled  Data  Variable  Structure  Control .  3 

2.2  Linear  Systems .  5 

2.2.1  Exogenous  Disturbances .  7 

2.2.2  System  Parameter  Variations  .  8 

2.2.3  Control  Coefficient  Matrix  Variations .  8 

2.3  Nonlinear  Systems  .  10 

2.3.1  Matched  Disturbances .  11 

2.3.2  Control  Coefficient  Variations .  12 

2.4  Stochastic  Systems .  14 

2.4.1  Continuous  Measurement .  15 

2.4.2  Discrete  Measurement .  17 

2.5  Experimental  Example  -  0{T^)  Sliding  Mode  Control  on  a  Flexible 

Structure  .  18 

3  VIBRATION  SUPPRESSION  IN  FLEXIBLE  STRUCTURES  VIA 

GENERALIZED  SLIDING  MODES  23 

3.1  Generalization  of  the  Sliding  Mode  (.'ontrol  Concept .  24 

3.2  Block  Control  of  Differential-Difference  Systems .  25 

3.3  Sliding-Mode  Control  for  Suppressing  Vil)rations  of  a  Flexible  Bar  .  .  27 

3.4  Dispersive  Flexible  Structures:  an  Integral  Transformation  Approach  30 

3.5  Use  of  Neural  Networks  with  the  Integral  Transformation  Approach  .  31 

4  RADIAL  BASIS  FUNCTION  NEURAL  NETWORKS  AND  WAVELET 

NEURAL  NETWORKS  35 

4.1  Radial  Basis  Function  Networks .  35 


4.2  Wavelet  Neural  Network  Architecture  1 .  37 

4.3  Multiresolution  Signal  Analysis .  38 

4.4  Wavelet  Neural  Network  Architecture  2 .  40 

4.5  Wavelet  Neural  Networks .  41 

4.6  Network  Design .  42 

4.7  RBF  Neural  Network  Simulation  Results . . .  46 

4.8  Wavelet  Neural  Network  Simulation  Results .  47 

5  MODELING  AND  CONTROL  OF  STRUCTURAL  SYSTEMS  US¬ 
ING  WAVELETS  53 

5.1  Background  on  Time-Frequency  Wavelets .  54 

5.2  Structural  System  Modeling  and  Identification .  59 

5.3  Multiresolutional  Sensitivity  Minimization .  61 

5.4  Multiresolutional  Controller  Development  and  Implementation  ....  63 


6  CONCLUSION 


71 


List  of  Figures 


2.1  Discrete  time  sliding  mode  in  sampled-data  systems .  4 

2.2  LIVE  truss  structure .  19 

2.3  Frequency  content:  Sliding  mode  control  (c=:l)  v.s.  No  control  ....  21 

2.4  Frequency  content:  Sliding  mode  control  (c=10)  v.s.  No  control  .  .  .  22 

3.1  Configuration  for  Training  Neural  Network .  32 

3.2  Closed  Loop  System  with  Neural  Network .  33 

4.1  RBF  and  Wavelet  NN  Architecture  for  Learning .  42 

4.2  Grid  of  wavelets  at  multiple  resolutions .  43 

4.3  Scatter  Plot  of  Nearest  Neighbor  Distance  for  NN  Training  Data  ...  46 

4.4  RBFNN  De.sign  4  -  Testing  Run  1 .  49 

4.5  RBFNN  Design  4  -  Testing  Run  2 .  50 

4.6  RBFNN  Design  4  -  Testing  Run  4 .  50 

4.7  WNN  Design  6  -  Testing  Run  1 .  51 

4.8  WNN  Design  6  -  Testing  Run  2 .  51 

4.9  WNN  Design  6  -  Testing  Run  4 .  52 

5.1  Partition  of  intervals,  {h},  constructed  using  disjoint  action  regions 

at  the  interval  endpoints .  57 

5.2  Typical  Time- Frequency  Wavelet .  58 

5.3  Example  of  possible  splitting  locations  for  a  particular  signal  and  a 

potential  orthogonal  basis  and  a  potential  level  basis  .  65 

5.4  Procedures  for, performing  cosine  packet  analysis  and  synthesis  on  a 

discrete  signal,  x{k),  for  a  particular  depth  level,  d .  66 

5.5  Impulse  and  frequency  response  for  a  flexible  structure  with  two  modes  66 

5.6  Cosine  packet  decomposition  and  best  luisis  selection  for  a  flexible 

structure  with  two  modes  at  cui  =  27r  and  uj2  =  47r  with  the  identi¬ 
cal  damping  factor  of  ^  =  0.01 .  67 

5.7  Reconstruction  of  the  impulse  response  for  a  flexible  structure  with 

two  modes .  68 

5.8  Impulse  and  frequency  response  for  a  flexible  structure  with  two  modes  68 


5.9  Cosine  packet  decomposition  and  best  level  basis  selection  for  a  flexible 

structure  with  two  modes  at  wi  =  2it  and  u:2  =  47r  with  the  identical 
damping  factor  of  C  =  0.01 . 

5.10  Reconstruction  of  the  impulse  response  for  a  flexible  structure  with 

two  modes . 

5.11  Feedback  Control  System . 


X 


List  of  Tables 


4.1  RBFNN  Training  Design  Summary .  47 

4.2  RBFNN  Testing  Design  Summary .  47 

4. -3  Grid  of  Wavelets  for  Actuator  Data .  48 

4.4  WNN  Training  Design  Summary  .  48 

4.5  WNN  Testing  Design  Summary .  49 


XI 


1.  INTRODUCTION 


In  this  report,  we  describe  the  total  effort  at  The  Ohio  State  University  on  the  Project 
Analysis  and  Control  of  Interconnected  Structures,  initiated  on  September  1, 
1992,  sponsored  by  the  Air  Force  Office  of  Scientific  Research  (AFSC),  under  Contract 
F49620-92-J-0460. 

The  main  thrust  of  our  effort  has  been  to  consider  the  control  of  interconnected  rigid 
and  flexible  structures.  In  order  to  accomplish  this,  we  have  utilized  different  analysis 
approaches  including  partial  differential  equations  representing  simple  substructures 
and  finite  element  models,  radial  basis  function  neural  networks  and  wavelet  neural 
networks  representing  complex  structures.  Finally,  we  have  introduced  the  use  of  a 
multiresolutional  technique  with  wavelets.  .4s  far  as  control  techniques  are  concerned, 
we  have  concentrated  on  sliding  mode  control  with  a  detailed  analysis  of  sampling 
effects.  The  multiresolutional  wavelet  approach  which  we  concentrated  on  during  the 
last  year  has  naturally  led  to  a  multiresolutional  controller  design  in  the  domain. 

The  following  chapters  will  provide  details  of  the  results  obtained  through  the  tenure 
of  the  project. 

In  Chapter  2  several  results  for  sliding  mode  control  of  sampled-data  systems  will 
be  outlined.  The  sampled-data  sliding  mode  control  results  will  be  used  to  design 
controllers  for  a  flexible  structure.  These  controllers  were  verified  on  an  actual  flexible 
structure  (which  was  developed  under  a  previous  AFOSR  grant). 

In  Chapter  3,  we  introduce  a  new  approach  for  vibration  damping  in  flexible  struc¬ 
tures  based  on  the  sliding-mode  control  approach.  Previous  use  of  this  approach 
has  been  for  flnite  dimensional  (approximate)  models  of  flexible  structures.  Here,  we 
retain  the  infinite-dimensional  model  for  the  structures  and  investigate  exact  solu¬ 
tions  using  sliding-mode  control.  The  use  of  neural  networks  in  conjunction  with  this 
method  are  also  investigated. 

In  Chapter  4,  consideration  of  neural  networks  let  to  our  development  of  a  new 
general  design  procedure  for  radial  basis  function  neural  networks  and  wavelet  neural 
networks.  These  type  of  neural  networks  are  especially  well-suited  for  dealing  with 
irregularly  positioned  data  points.  We  present  simulation  results  for  these  two  types  of 
networks  based  on  predicting  actuator  nonlinearities  occurring  in  a  flexible  structure. 

In  Chapter  5,  we  explore  uses  for  time-frequency  wavelets.  These  include  the  mod¬ 
eling  and  control  of  structural  systems  and  subsystems.  We  develop  a  versatile  struc¬ 
ture  model  based  on  time-frequency  wavelets  as  opposed  to  the  traditional  finite 
element /modal  models.  The  wavelet  models  are  developed  for  the  purpose  of  being 
utilized  during  the  controller  design  phase.  We  then  present  some  theoretical  results 
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necessary  for  developing  a  sensitivity  minimization  controller  based  on  a  wavelet 
model. 

This  project  is  also  linked  to  AASERT  Project  F49620-92-J-0299  entitled  Integrated 
Circuits  for  Distributed  Control-(AASERT  FY91).  The  main  results  of  the  research 
performed  by  the  students  supported  under  this  AASERT  Project  are  documented 
separately,  although  some  of  the  neural  network  related  efforts  are  reported  here  for 
completeness.  (The  AASERT  document  concentrates  on  related  yet  distinct  efforts  of 
developing  and  using  a  FEM  analog  integrated  circuit  chip  for  substructure  control.) 


2.  IMPLEMENTATION  OF  VARIABLE  STRUCTURE  CONTROL 
FOR  SAMPLED-DATA  SYSTEMS 


2.1  Background  on  Sampled  Data  Variable  Structure  Control 

The  characteristic  feature  of  a  continuous  time  variable  structure  systems  (VSS)  is 
that  sliding  mode  occurs  on  a  prescribed  manifold,  or  switching  surface,  where  switch¬ 
ing  control  is  employed  to  maintain  the  state  on  that  surface  [1,  2.  3,  4,  5].  Since  the 
theory  has  been  originally  developed  from  a  continuous  time  prospective,  implemen¬ 
tation  of  sliding  mode  for  sampled-data  systems  encounters  several  incompatibilities 
due  to  limited  sampling  rate,  sample/hold  effect,  and  discretization  errors.  As  a 
result,  a  direct  translation  of  continuous  time  variable  structure  control  design  for 
discrete  implementation  leads  to  the  chattering  phenomenon  in  the  vicinity  of  the 
switching  surface. 

This  chapter  deals  with  the  implementation  of  variable  structure  control  for  sampled- 
data  systems  by  maintaining  sliding  mode  in  discrete  time.  Although  a  considerable 
amount  of  work  has  been  done  analyzing  discrete  time  sliding  mode,  very  little  was 
directly  addressed  to  the  sampling  issues.  In  [6],  Milosavljevic  studied  the  oscillatory 
characteristic  (quasisliding)  in  the  neighborhood  of  the  discontinuity  planes  due  to 
discretization  of  continuous  time  signals.  Existence  conditions  of  quasisliding  mode 
were  derived  as  a  discrete  extension  from  the  continuous  time  VSS  theory.  For  con¬ 
trol  law  design,  Utkin  and  Drakunov  proposed  a  definition  of  discrete  time  equivalent 
control  that  directs  the  states  onto  the  sliding  surface  in  one  sampling  period  [7],  To 
remain  on  the  surface,  the  associated  control  appears  to  be  non-switching.  Subse- 
cpiently,  the  theoretical  basis  was  furnished  with  a  formal  definition  of  sliding  mode 
for  discrete  time  systems  in  the  context  of  semigroups  [4]. 

Related  works  from  different  prospectives  can  be  found  in  [8,  9,  10,  11,  12,  13].  Sarp- 
turk  et  aL  took  a  Lyapunov  point  of  view  for  discrete  time  linear  systems.  It  was 
asserted  that  the  switching  control  be  bounded  in  an  open  interval  to  guarantee  con¬ 
vergence  of  sliding  motion.  This  interval  was  later  found  to  depend  linearly  on  the 
distance  of  the  state  from  the  switching  surface  [9],  which  suggested  a  nonswitch¬ 
ing  control  when  discrete  time  sliding  mode  is  attained.  Sira-Ramirez  imposed  the 
geometrical  concept  to  general  nonlinear  SISO  systems  along  with  the  existence  is¬ 
sues  [11].  Parallel  results  with  continuous  time  VSS’s  were  obtained.  To  implement 
discrete  time  sliding  mode  control  law  in  compliance  with  the  existence  conditions, 
Furuta  proposed  the  idea  of  sliding  cone  (sliding  sector),  where  switching  control 
takes  place  only  when  the  states  are  out  of  the  sector,  while  inside  of  the  sector, 
the  control  law  remains  continuous  [10].  Kaynak  and  Denker  used  an  ARM  A  model 
to  characterize  the  control-sliding  surface  relationship  and  yielded  a  non-switching 
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Sliding  surface  s(x)=0.  _ _ ^ _ _ _ _ _ 
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Figure  2.1:  Discrete  time  sliding  mode  in  sampled-data  systems 

type  of  control  with  a  predictive-corrective  scheme  [12].  Bartolini  et  al.  incorporated 
adaptive  control  strategy  to  account  for  system  uncertainties  and  designed  control 
law  in  terms  of  a  discrete  time  equivalent  control  [13]. 

Our  work  has  mainly  extended  the  original  idea  of  [7]  to  more  general  problems  with 
detailed  analysis  of  the  sample/hold  effects.  Preliminary  results  were  reported  in  [14, 
15,  16]  and  the  full  development  is  presented  here.  The  control  objective  is  to  maintain 
the  states  on  the  switching  surface  at  each  sampling  instant.  Between  samples,  the 
states  are  allowed  to  deviate  from  the  surface  instead  of  being  constantly  and  exactly 
on  the  switching  surface,  even  the  equivalent  state  travels  within  a  boundary  layer 
of  that  surface  (see  Figure  2.1).  It  will  be  demonstrated  that  the  thickness  of  this 
boundary  layer  can  be  reduced  significantly  by  proper  consideration  of  the  sampling 
phenomenon  in  the  control  design.  Three  classes  of  systems  (linear,  nonlinear,  and 
stochastic  systems)  will  be  investigated  individually.  Robustness  against  both  internal 
and  external  uncertainties  will  be  considered.  Furthermore,  the  chattering  problem 
will  be  addressed. 

The  chapter  is  organized  as  follows:  in  Section  2.2,  we  consider  discrete  time  slid¬ 
ing  mode  in  linear  sampled-data  systems.  The  iritersample  behavior  as  well  as  the 
matching  condition  will  be  examined.  Robustness  against  system  uncertainties  are 
elaborated.  Three  types  of  uncertainties  will  be  studied  individually.  The  overall 
effects  can  be  treated  as  a  linear  combination  of  the  three.  In  Section  2.3,  the  re¬ 
sults  are  extended  to  nonlinear  systems.  Based  on  similar  assumptions,  the  control 
strategies  can  be  directly  translated  from  the  linear  case.  Section  2.4  deals  with  linear 
stochastic  systems.  Sliding  mode  is  defined  in  terms  of  a  conditional  probability  func¬ 
tion,  given  the  cr-algebra  generated  by  the  measured  state  information.  Both  cases 
with  continuous  measurement  and  discrete  measurement  will  be  considered.  To  min¬ 
imize  the  deviation  from  the  sliding  manifold  in  the  mean  square  sense,  the  optimal 
filtering  problems  will  be  solved.  In  each  of  tlu'  above,  the  chattering  phenomenon 
will  be  effectively  removed  with  robustness.  .\u  ('xperimental  example  of  a  flexible 
structure  vibration  control  is  given  in  section  2.5  to  demonstrate  the  effectiveness  ot 
the  proposed  sampled-data  sliding  mode  control  technique.  Finally,  the  conclusion 
section  summarizes  the  chapter. 
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2.2  Linear  Systems 


Consider  the  dynamic  system  with  a  prescribed  switching  surface 

X  =  Ax  +  Bu  +  Df 
<5  =  {a;|s(ar)  =  Cx  =  0}, 

where  the  state  x  G  the  control  u  G  jR’”,  the  disturbance  vector  /  G  and 
the  sliding  surface  vector  s  G  i?™;  A,  B,  D  are  constant  matrices  of  appropriate 
dimensions;  and  the  m  x  n  matrix  C  is  chosen  so  that  the  system,  when  traveling 
on  S,  will  achieve  the  desired  sliding  mode  dynamics  [2],  VVe  say  the  disturbance  is 
“rejectable”  if  there  exists  a  control  u  that  can  eliminate  its  effect  instantaneously  in 
the  dynamic  system.  If  the  closed  loop  system  behavior  does  not  depend  on  /(t),  we 
say  the  controlled  system  is  disturbance  invariant.  It  is  well  known  that  the  matching 
condition  [3] 

rank[B,  D]  =  rank[B]  (2.3) 

has  to  be  satisfied  for  the  system  ((2.1))  to  be  disturbance  invariant  in  sliding  mode. 

The  discrete  time  representation  of  the  dynamic  equation  ((2.1))  is  obtained  by  ap¬ 
plying  u  through  a  zero-order- hold: 

.Tfc+i  =  ^Xk  +  Tuk  +  dk,  (2.4) 

where  $  =  F  =  /J"  e^^dXB,  dk  =  Jq  e‘‘^^Df{(k  +  1)T  —  \)dX.  Note  that  the 
magnitude  of  F  and  dk  are  both  0{T)  if  f{t)  is  bounded. 

The  discrete  time  sliding  mode  control  is  to  steer  the  states  towards  and  maintain 
them  on  the  surface  S  at  each  sampling  instant  such  that 

Sk  =  Cxk  =  0  k  >  ks.  (2.5) 

During  the  sampling  interval  kT  <  t  <  {k  +  [)T.  the  state  will  deviate  from  S.  The 
following  Lemma  provides  bounds  on  the  deviation  from  the  sliding  surface. 

Lemma  1  If  the  disturbance  f{t)  in  ((2.1))  is  Ijounded  and  smooth  in  the  intersample 
interval  kT  <  t  <  {k  +  l)T,  then  <  ^  implies  s{t)  =  0{T^)  during  that 

(  s/j+i  —  U 

sampling  period. 
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The  0{T^)  deviation  is  inevitable  for  sampled-data  systems.  It  is  the  ultimate  per¬ 
formance  a  sampled  controller  can  attain  for  continuous  time  plants.  Note  that  if 
Sk  —  i'T^)  and  =  0{T^),  the  intersample  value  of  the  sliding  surface  vector  is 
still  0{T^). 

The  sample/hold  effect  not  only  destroys  perfect  sliding  mode  as  depicted  above, 
it  also  invalidates  the  complete  disturbance  rejection  property  for  the  discrete  time 
system  ((2.4)). 

Lemma  2  The  sampled-data  system  ((2.4))  is  disturbance  invariant  on  the  discrete 
time  sliding  manifold  ((2.5))  only  if  the  continuous  time  matching  condition  ((2.3)) 
holds. 

The  above  lemma  states  the  necessity  of  the  continuous  time  matching  condition, 
however,  the  sufficiency  is  not  confirmed.  The  zero-order-hold  is  applied  to  the  control 
variables  only.  A  similar  “hold”  does  not  take  place  in  the  disturbance  channels.  In 
general,  the  discrete  time  disturbance  dk  will  not  lie  in  the  range  space  of  the  control 
coefficient  matrix  T  and  hence  can  not  be  rejected  completely  even  if  discrete  time 
sliding  mode  occurs. 

Although  perfect  sliding  mode  and  complete  disturbance  rejection  are  not  possible  in 
sample-data  systems,  one  can  still  maintain  the  states  in  the  vicinity  of  the  sliding 
surface  and  retain  satisfying  disturbance  rejection  character  by  proper  consideration 
of  the  sampling  phenomenon  in  the  control  design. 

We  will  consider  three  different  types  of  uncertainties  that  are  important  in  affecting 
a  dynamic  system’s  behavior;  namely,  exogenous  disturbances,  system  parameter 
variations,  and  control  coefficient  variations.  It  will  be  shown  that  these  uncertainties 
can  be  rejected  to  at  least  0{T^)  accuracy  and  that  the  chattering  phenomenon 
inhabiting  in  continuous  time  variable  structure  systems  can  be  removed. 

Since  the  control  in  the  sampled  data  system  will  show  jumps  from  the  continuous 
time  point  of  view,  the  formal  concept  of  continuity  is  not  compatible  with  the  notion 
of  discrete  time  control!  In  order  to  make  the  analogy  in  between,  let  us  inti'oduce 
the  following  definition: 

Definition  1  The  discrete  time  control  law  Uk  is  said  to  be  ecjuivalent  to  discontin¬ 
uous  if  Aufc  =  0(1),  continuous  if  Aufc  =  0(7),  smooth  if  A^n*,  =  0(7^),  where  A 
denotes  the  backward  differencing  operator  1  — 

To  avoid  a  cumbersome  description,  the  “continuity”  of  discrete  time  control  to  appear 
in  the  following  text  will  refer  to  its  equivalent  meaning  as  defined  above. 


6 


2.2.1  Exogenous  Disturbances 


Consider  the  case  when  only  the  exogenous  disturbance  f{t)  is  present.  We  define 
the  discrete  time  equivalent  control  by  solving  Sk+i  =  0  [7],  which  leads  to 

ul^  =  -{CT)-^C{^Xk  +  dk)  (2.6) 

assuming  CF  is  invertible.  Note  that  ul?  is  generally  not  accessible  because  it  requires 
the  value  of  dk,  which  depends  on  the  future  values  of  the  disturbance  f{t)  from  the 
present  sampling  instant  kT  to  the  next  one  {k  +  l)r.  It  is  impossible  to  evaluate 
dk  exactly  unless  f{t)  is  known.  Nevertheless,  with  boundedness  and  smoothness 
assumptions  imposed  on  /(t),  dk  can  be  predicted  by  its  previous  value  dk-i  [17], 
which  can  be  computed  from  ((2.4)).  The  error  is 

fT  /-(fc+ijr-A  . 

dk-dk-i=  f(a)dadX  =  0{T^) 

Jo  J  kT — A 


assuming  f{t)  is  bounded.  To  approximate  we  define 

ur  =  -(CT)-iC($Xfc  +  4-i).  (2.7) 

Taking  the  Taylor’s  series  expansion  of  $  yields  the  relationship  C^Xk  =  5^  +  [AT  + 
•  •  •)xk.  Since  the  magnitude  of  CT  and  C4— i  are  both  0{T),  in  order  that  i4" 
be  admissible,  the  magnitude  of  C^Xk  must  also  be  0[T).  In  other  words,  Xk  must 
lie  within  an  0{T)  boundary  layer  of  S  (or  Sk  =  0{T))  for  n)”  to  be  admissible. 

In  order  to  define  the  boundary  layer  of  the  sliding  surface,  let  us  adopt  the  following 
notation: 

St  =  {xkWk  €  7/}  is  the  0[T)  boundary  layer  of  S. 

We  propose  the  variable  structure  control  law 


—  (CT)  ^C{^Xk-\- dk-i)  when  Cj, 
Uk  =  <  —{CV)~^C[^Xk  +  dk-i  —  Sk  +  A  ■■<(jn(sk))., 

otherwise. 


(2.8) 


The  positive  definite  matrix  K  will  determine  the  stepsize  for  the  state  to  approach 
the  boundary  layer  Sj-  The  magnitude  of  K  has  to  be  chosen  small  enough  not  to 
overshoot  St- 
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Theorem  1  For  a  linear  samplecl-data  system  ((2.4))  with  the  exogenous  disturbance 
satisfying  the  matching  condition  ((2.3)),  the  effect  of  f(t)  can  be  reduced  to  0{T'^) 
if  the  variable  structure  control  law  ((2.8))  is  applied. 

2.2.2  System  Parameter  Variations 

Consider  a  dynamic  system  subject  to  parameter  variations 

i  =  {A  + /\A{t))x  +  Bu.  (2.9) 

The  matching  condition  for  the  uncertainty  A/l(f)  is 

Rank[AA{t)Q,  B]  —  Rank[B]  Vt,  (2.10) 

where  0  is  an  n  x  (n  —  m)  matrix  with  the  columns  forming  a  basis  of  the  subspace 
nuU(C)  [3].  The  unknown  AA(t)a:  may  be  treated  as  a  disturbance.  The  associated 
sampled-data  representation  is  given  in  ((2.4)),  where 

4  =  r  e'^^AAiik  +  1)T  -  X)x{{k  +  1)T  -  X)dX.  (2.11) 

Jo 

Assuming  AA{t)  is  bounded,  cR  can  be  predicted  by  4-1  with  O(T^)  error,  and  we 
have  the  following  theorem. 

Theorem  2  For  a  linear  system  with  parameter  variations  ((2.9)),  the  control  ((2.8)) 
will  lead  to  s^+i  =  0(T^)  if  AA{t)  satisfies  the  matching  condition  ((2.10)).  Further¬ 
more,  the  state  will  attenuate  to  zero  asymptotically. 

2.2.3  Control  Coefficient  Matrix  Variations 

Systems  with  unknown'  control  coefficient  matrix  variations  have  the  following  dy¬ 
namic  equation 

i  =  Ax  +  {B  +  AB{t))u.  (2.12) 

The  matching  condition  is 

Rank[AB{t),  B]  =  Rank[B]  Vt.  (2.13) 

The  system  dynamic  equation  in  discrete  time  is  the  same  as  ((2.4)),  where 
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dk  =  ATkUk, 

AFfc  =  e^'^AB{ik  +  1)T  -  X)dX. 


Employing  the  control  law  u‘^  induces  a  dynamic  feedback  in  Uk  through  f4-i-  Since 

Sk+i  =  c^xk  +  C(r  +  Ark)uk,  {2-14) 

the  stability  of  Sk  implies  stability  in  Xk,  which  ensures  stability  in  Uk  if  ^(r  +  AP^) 
is  invertible. 

Introduce  the  following  notation: 

=  AVk  -  ATk-u  A^  =  C’Ar,(Cr)-E  (2.15) 

The  qiiantit}^  can  be  regarded  as  a  measure  of  variations  in  the  control  coefficient 
matrix  with  respect  to  the  known  one  F.  If  AB{t)  and  AB{t)  are  bounded,  we  have 
—  0{T^)  and  Ak  =  0(1)- 

Lemma  3  For  the  linear  system  with  control  coefficient  matrix  variation  ((2.12)),  if 
AB{t)  satisfies  the  matching  condition  ((2.13));  and  the  states  Xk-i^  Xk  are  in  the 
boundary  layer  <St,  then  the  control  t/JF  will  lead  to 

+  Ak-^iSk^x  +  (2.16) 


with  Sk  =  0{T'^). 

Stability  of  Sk  in  ((2.16))  can  be  assured  if  Ak  is  small  enough;  that  is,  the  variation 
of  AB{t)  is  small  with  respect  to  B.  Since  the  quantity  6k  is  O(T^),  the  state  will 
enter  the  0{T^)  boundary  layer  of  S  if  ((2.16))  is  stable. 

Theorem  3  The  stdite  xk  attenuates  to  zero  asymptotically  if 

L  The  conditions  in  Lemma  3  hold. 

2.  Equation  ((2A6))  is  stable, 

3.  The  m  X  m  matrix  C{T  +  ATk)  is  no  using  alar. 

If  all  three  types  of  uncertainties  are  present  in  the  dynamic  ecjuation,  the  variable 
structure  control  law  ((2.8))  can  still  be  employed  to  compensate  the  uncertainties 
simultaneously.  In  this  case,  dk  is  the  lumped  effect  of  all  disturbances,  the  resultant 
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sliding  mode  accuracy  is  s{t)  =  0{T^).  As  long  as  the  states  remain  in  the  0{T'^) 
boundary  layer  of  the  sliding  surface,  becomes  continuous  and  chattering  in  the 
control  is  removed. 

Theorem  4  The  control  is  continuous  if  =  0{T^).  Furthermore,  ui^  ap¬ 
proaches  the  continuous  time  equivalent  control  as  the  sampling  rate  approaches 
infinity. 

uk  ^  =  -{CB)-^C{Ax  +  Df).  (2.17) 

2.3  Nonlinear  Systems 

In  the  previous  sections,  we  have  developed  a  discrete  time  sliding  mode  control 
method  applied  to  linear  sampled-data  systems  to  alleviate  chattering.  Here,  we 
extend  it  to  nonlinear  systems.  We  assume  the  nonlinearities  and  uncertainties  are 
bounded  and  smooth,  in  this  case  the  past  values  of  the  states  can  be  utilized  to 
predict  the  effect.  This  approach  allows  us  to  keep  the  state  in  the  0{T^)  vicinity  of 
the  sliding  manifold. 

Consider  the  nonlinear  system 

X  —  a(x)  +  b{x)u  -f  v(x,t)  (2.18) 

with  a  prescribed  sliding  manifold 

S  =  {a:|s(;r)  =  0},  (2.19) 

where  a(x)  and  b(x)  are  known  differentiable  functions  with  dimensions  n  x  1  and 
n  X  m  respectively;  v{x,t)  is  an  n  x  1  unknown  vector-valued  differentiable  function. 
The  disturbance  matching  condition  is 

Rank[b{x),v{x,t)]  =  R(i,nk[b{x)]  when  .r  G  >S.  (2.20) 

Let  G{x)  =  ■§^s{x)  be  the  m  x  n  Jacobian  matrix  of  .s(,r)  so  that  the  dynamics  of  the 
sliding  surface  vector  can  be  expressed  as 

s{x)  =  G{x)a{x)  -f  G{x)b{x)u  +  G{x)v(x.t).  (2.21) 

Applying  a  zero-order-hold  to  the  control  yields  the  discrete  time  representation  of 
((2.18))  and  ((2.21)) 
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Xk+1  =  Xk  +  a{x)dt  +  b{x)dt  Uk 

+  v{x,t)dt 


sk+i  =  G{x)a{x)dt  +  G{x)b{x)dt  Uk  ,,  .,,, 

+  G{x)v{xG)dt. 

The  control  objective  is  to  keep  the  states  Xk  on  the  sliding  surface  <5  in  discrete  time 
s{xk)  =  0,  for  k  >  kg  (2.24) 

The  following  Lemma  is  the  nonlinear  version  of  Lemma  2. 

Lemma  4  The  sampled-data  system  ((2.22))  is  disturbance  invariant  on  the  discrete 
time  sliding  manifold  ((2.24))  only  if  the  continuous  time  matching  condition  ((2.20)) 
holds. 


2.3.1  Matched  Disturbances 

Consider  the  system  ((2.18)),  in  which  the  disturbance  v{x,t)  satisfies  the  matching 
condition  ((2.20)).  Introduce  the  notation: 

Fk  =  G{x)a{x)dt 

Bk  =  G{x)bix)dt 

hk  =  G{x)v{xG)dt 


SO  that  the  discrete  representation  ((2.23))  can  be  rewritten  as 


^  Fk  +  Bk  Uk  hk<  (2.26) 

The  equivalent  control  is 

u^k  =  -FZ^i^k  +  Fk  +  hk)  (2.27) 

assuming  Bk  is  invertible.  The  magnitude  of  Fk-,  Bk,  and  hk  is  0{T)  if  G{x),  a{x). 
b{x),  and  v{x,t)  are  bounded.  Fk,  Bk,  and  hk  usually  cannot  be  evaluated  exactly 
since  they  depend  on  the  future  value  of  the  state  and  the  unknown  disturbance. 
Nevertheless,  with  some  mild  smoothness  condition  on  the  functions  G{x),a{x),b{x) 
and  u(x,f),  they  can  be  approximated  by  the  following  expressions: 
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(2.28) 


Fk  =  G{xk)a{xk)T 
Bk  =  G{xk)bixk)T_ 
hk  =  Sk+i  —  Sk  —  Fk  —  Bk. 


With  ((2.28)),  we  obtain  an  approximation  for  the  equivalent  control  ((2.27)) 

^  +  hk-\)^  (2.29) 

where  is  invertible  by  assumption.  Similar  to  ((2.7)),  must  lie  in  the  boundary 
layer  St  for  to  be  admissible.  The  proposed  variable  structure  control  law  is: 

„  _  (  -Bj;\sk  + Fk  +  Jik-i)  _  when  X‘^^e  St 

^  \  -B^^{Ksgn{sk)  +  Fk  +  hk-i)  otherwise 

Theorem  5  If  the  disturbance  in  dynamic  system  ((2.18))  satisfies  the  matching 
condition  ((2.20)),  then  the  effect  of  v{x,t)  can  be  reduced  to  0{T)  if  the  control 
((2.30))  is  applied. 


2.3.2  Control  Coefficient  Variations 

Consider  the  case  when  there  is  an  unknown  variation  b{x)  in  the  control  coefficient 
X  =  a(x)  +  b(x)u  +  b{x)u.  (2.31) 

The  matching  condition  is 

Rank[b{x),  6(a;)]  =  B.ank[b{x)]  when  x  €  S.  (2.32) 

Following  the  same  discretization  procedure  yields  the  discrete  representation  indi¬ 
cated  in  ((2.26))  with 

hk  =  BkUk 

Bk  =  Gix)b{x)dt. 

Introduce  the  following  notation: 

ep^  =  Fk  -  Fk,  ep^  -  ~  -  Fk-\, 

=  Bk  -  Bk,  -  Bk  -  Bk-\, 

eBk  =  Bk  -  Bk-\,  ep^  =  Bk  -  Bk-\,  (2.33) 

c/i;.  —  hk  hk—i,  hk  hk 

Ak  =  BkB^^ ,  =  Afc  -  Ak-\. 
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It  can  be  proved  that  A*,  =  0(1),  =  0(T),  and  ej^.  =  0{T'^),  =  O(T^),  where 

j  =  F,  B,  F,  B,  B. 

Lemma  5  When  the  state  is  in  the  boundary  layer  St,  the  following  equalities  are 
true: 


•5fc+i  —  hk  —  hk-i  +  al{2sk  —  Sk-\]  +  ctl  (^-34) 

hk  =  —  AA;[(sfc  +  Fk  F  hk-i)  +  /3l{2sk  —  Sk-i)  +  (2.35) 

where  aj  =  0{T),  al  =  0(T^},  /JJ  =  0(7),  =  0{T^). 

The  equalities  in  Lemma  (5)  can  be  obtained  by  direct  substitution  from  ((2.33))  and 
((2.29))  into  ((2.26)).  Furthermore, 

al  =  -2es^By 

=  -^Bk^k  ^Fk  +  ep^  +  epk^k 

F^B^^k  +  Bk-iUk-i)  (2.36) 

/3l  =  eB,B;^ 

l^k  =  ~^Fk  ~  ^hk-i^Bk^k  ~  Bk^iUk-i). 

They  serve  as  preparation  for  the  theorem  stated  below: 

Theorem  6  For  the  nonlinear  system  with  control  coefficient  matrix  variation  ((2.31)), 
the  control  ((2.30))  will  lead  to  Sk  — >  0{T^)  asymptotically  if 

1.  The  state  lie  in  the  boundary  layer  Sj- 

2.  The  variation  b{x)  satisfies  the  matching  condition  ((2.32))  and  that  b{x)  is 
bounded. 

3.  The  following  dynamic  equation  for  Sk  is  stable. 

•sa;+i  =  +  a^Sfc_i  +  a\sk-2  +  bk,  (2.37) 


where 
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4  =  -2Afc  -  2A,/?,i  +  2al  =  -2A,  +  0(T) 

4  =  Ak-i  +  2Akal  —  +  Ak,dl:  +  'lAk-i^  —  a[ 

=  A,_x  +  0{T) 

al  =  -  A, 4  -  A,_i4  =  0(T) 

h  =  —  A^.e^fc  —  €i\^Fk-i  +  Akol 
+  ^Aki-^f:-2  +  Bk-2Uk-2) 

-Ak4  +  ^k-ifiU  +  4  =  0[T^). 

If  the  sampling  rate  is  sufficiently  high,  the  dynamic  equation  ((2.37))  is  reduced  to 
second  order 


■sjt+i  —  —‘2AkSk  +  Ak-iSk-i-  (2.38) 

This  agrees  with  ((2.16)),  the  result  for  linear  systems.  Stability  of  Sk  can  be  assured 
if  Ak  is  small  enough;  that  is,  the  variation  of  the  control  coefficient  g{x)  is  small 
with  respect  to  b{x).  If  the  0{T)  terms  in  ((2.37))  are  not  negligible,  the  dynamic 
equation  for  Sk  becomes  third  order  with  an  0{T^)  input.  Stability  of  ((2.37))  implies 
that  Sk  — >  0{T^)  asymptotically. 

2.4  Stochastic  Systems 

This  section  is  devoted  to  discrete  time  sliding  mode  control  for  continuous  time  linear 
systems  with  stochastic  disturbances.  The  disturbance  is  assumed  to  satisfy  an  Ito 
type  stochastic  differential  equation  [18] 

dx  =  Axdt -\- Budt -\- [D^dt -\- Rdw\)  (2.39) 

di  =  Widt  +  QdiV2,  (2.40) 

where  tr>i(t),  W2{t)  are  independent  standard  Wiener  processes.  Here  represents 
the  colored  part  of  the  noise  and  the  formal  expression  the  white  part.  Without 

loss  of  generality,  the  autocovariance  matrices  of  t  he  white  noises  are  assumed  to  be 
identity. 

We  consider  two  possible  state  information  channels;  continuous  measurement  and 
discrete  measurement.  To  compensate  for  the  stocliastic  disturbances  with  a  sampling 
controller,  a  mixed  continuous-discrete  type  of  liltering  problem  will  be  solved  it 
continuous  state  information  x[t)  is  used.  On  I  lie  other  hand,  the  use  of  discrete 
state  information  Xk  will  lead  to  a  discrete  linu'  optimal  filtering  problem.  In  both 
cases,  we  apply  the  following  assumption: 


Assumption  1 


1.  W  has  asymptotically  stable  eigenvalues; 

2.  R  is  full  rank; 

3.  The  processes  Rwi{t)  and  Qw2{t)  are  uncorrelated. 

2.4.1  Continuous  Measurement 

Consider  the  discrete  time  version  of  ((2.39)) 

Xk+i  =  +  Tuk  +  i]k  +  Vk,  (2.41) 

where  r/k  =  is  colored  part  of  the  disturbance  and 

Vk  =  is  a  white  noise.  Since  the  white  noise  can  not  be 

compensated,  the  control  objective  is  to  minimize  the  effect  of  the  colored  part  pk  in 
the  sampled-data  system  ((2.41)).  Applying  the  control  law  with  an  estimation  of  pk 

uk  =  -{CT)-^c^xk  -  {cry^cfik,  (2.42) 

we  obtain 


Sk+i  =  C{pk  +  Vk  -  Pk)-  (2.43) 

For  Uk  ((2.42))  to  be  admissible,  it  still  requires  that  the  state  lie  in  the  0{T)  boundary 
layer  St-  This  boundary  layer  can  be  defined  similarly  as  in  the  previous  sections. 
Here,  we  only  consider  when  Xk  S  St- 

In  the  presence  of  stochastic  disturbances,  sliding  mode  is  defined  in  terms  of  a 
probability  measure.  To  minimize  £'||sfc+i||^  given  the  observation  of  x{t)  up  to  the 
instant  kT,  the  estimate  rjk  should  be  picked  as  a  conditional  expectation 

Pk  =  E{pk\TkT),  (2.44) 

where  J-t  denotes  a  cr-algebra  generated  l)y  the  process  x(t) 

Tt  =  (j{x(t)|0  <  t  <  t}.  (2.4.5) 

Since  wi{t)  is  a  martingale  [18],  the  estimation  i/k  is  obtained  by  taking  a  conditional 
expectation  given  the  cr-algebra  J-kT 
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f,k  =  /  =  Eiu, 

JkT 


(2.46) 


where  the  matrix  E  is  of  the  form 

=  =  de'^’-DdX, 

Jo 

and  if  is  a  conditional  expectation 

m  =  £(«i)|jf,). 


(2.47) 


(2.4S) 


Therefore,  from  ((2.46))  it  suggests  that  in  order  to  obtain  the  values  of  7]k  we  need 
a  filter  to  calculate  The  following  theorem  provides  the  form  of  such  a  filter. 

Theorem  7  If  assumption  (1)  holds,  the  optimal  estimate  of  ft  is 

ik  =  Zk  +  LR^Xk,  (2.49) 


where  Zk  satisfies 


Zk  —  ^Zk-i  +  (^k-i  +  nu^_i 


(2.50) 


with  A,  n  and  (k  defined  by: 

A  = 

a  =  ffe''UXM,  (2.51) 

=  /(f  ,)T<!"'*4'-‘IWi(A)dA, 

where  the  matrices  N,  M,  L  are 


H  =  W-  LEd'D, 

N  =  -LR^A  +  WLR'^  -  LR^DLR^, 
M  =  -LR?'B, 

L  =  PD^RiR'^R)-^ 


(2.52) 


and  P  is  the  positive  definite  solution  of  the  algebraic 


matrix  Ricatti  equation 


PW^  +  WP  -  PD^R{R^R)-^R^DP  +  QQ^  =  0. 


(2.53) 
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The  filter  ((2.49))  and  ((2.50))  provides  a  way  to  calculate  the  conditional  expection 
of  (k  given  the  information  from  the  continuous  time  process  x{t).  Substituting  the 
estimate  into  the  equation  ((2.46))  we  obtain  an  optimal  rjk  as: 

Vk  =  ’^{zk  +  LR^Xk).  (2.54) 

The  quantity  rjk  +  Vk  —  rjk  from  the  right  hand  side  of  ((2.43))  is  a  discrete-time  white 
noise  process  with  variance  of  order  0{T).  It  means  that  the  mean  sciuare  deviation 
of  from  zero  is  of  order  0(T^). 

2.4.2  Discrete  Measurement 

Suppose  we  only  have  access  to  the  discrete  state  information,  the  cr-algebra  generated 
by  the  discrete  measurement  Xk-,  Xk-\,  •  •  •  is 

:Fk  =  <  j  <  k).  (2.55) 

The  discrete  time  representations  of  ((2.39))  and  (??)  to  be  considered  are 

Xk+i  =  ^Xk  +  T^k  +  Tuk  +  Vk  +  Uk  (2.56) 

(2.57) 

where  i/fc  =  Jq  Jq  e^^^~^^Qdw2{kT+r)dX,  qk  = 

are  discrete  white  noises  generated  form  the  process  dw2\  the  matrices  'I'  and  T  are 

-X  =  e^'^  r  e-A^De^'^dX  (2.58) 

Jo 

To  minimize  the  effect  caused  by  ^k,  we  apply  tlie  control  law 

uk  =  -{cry^c^xk  -  {cvy^crik,  (2.59) 

where  ^k  is  the  optimal  estimation  of  (k  given  Xk 

ik  =  EUklJ^k).  (2.60) 

From  ((2.56)),  the  discrete  measurements  Xk  and  Xk-\  yield  the  innovation  process 

ijk-i  =  Xk  -  ^Xk-i  -  Viik-i  =  T^k-i  +  Vk-i  -f  i^k-i-  (2.61) 
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Combining  ((2.61))  with  ((2.57)),  one  can  construct  a  stationary  Kalman  filter  with 
the  process  noise  qk-\  in  ((2.57))  correlated  with  the  measurement  noise  (t'fc_i  + 
Vk-\)  in  ((2.61))  through  dw2-  To  obtain  a  standard  Kalman  filter  formulation  with 
uncorrelated  noises,  we  rewrite  equation  ((2.40)) 


^k+i  =  ^6  +  hr  +  qk, 


(2.62) 


where  H  and  R  are  the  covariance  matrices 


H  -  E{qk{vk  +  i^k)*),  R  =  E{{vk  +  i^k){vk  +  I'kT), 


and  qk  are  defined  as 

'5  =  ^  -  qk  =  qk  -  HR-\vk  +  Vk)- 

With  ((2.62))  and  the  innovation  process  an  optimal  estimation  for  ^k  is  acquired 


ik  =  #4_i  +  ^^-'yfc-i  +  Z(y;t-i-Ta-i),  (2.63) 

L  =  .PT^(TPT^  +  ^)-\  (2.64) 

P  =  t(P-PT^(TPT^  +  P)-^TP)$'^  +  OC^,  (2.65) 

where  is  the  autocovariance  matrix  of  the  discrete  white  noise  qk  in  ((2.62)). 
Applying  the  control  law  ((2.59))  yields 

•Sfc+i  =  C{'t^k  E  Vk  Vk  —  T(f/j).  (2.66) 


The  quantity  T((ffc  —  (ft)  +  Ufc  +  i/fc  from  the  right  hand  side  of  ((2.66))  is  also  a  discrete 
time  white  noise  with  variance  of  order  0{T).  Therefore,  the  mean  square  deviation 
of  Sk+i  from  the  sliding  surface  is  0{Ti). 

2,5  Experimental  Example  -  0{T^)  Sliding  Mode  Control  on  a  Flexible 
Structure 

This  section  deals  with  sampled-data  implementation  of  sliding  mode  for  a  flexible 
structure  control  problem.  We  apply  the  sampled-data  sliding  mode  control  tech¬ 
niques  to  a  flexible  structure  developed  for  Large  Interconnected  Vibration  Experi¬ 
ment  (LIVE)  at  the  Ohio  State  University.  The  truss  structure,  constructed  of  hol¬ 
low  PVC  tubes  and  Noryl  nodes,  is  aligned  on  a  vertically  cantilevered  configuration 
shown  in  Figure  2.2. 


IS 


Figure  2.2:  LIVE  truss  structure 


The  physical  properties,  control  hardware  setups,  preliminary  dynamic  modelling  and 
experimental  results  were  published  in  [19]  by  Redmill  and  Ozgiiner.  Extended  works 
in  sampled-data  sliding  mode  control  experiment  are  presented  here. 

The  dynamic  equations  of  the  vertically  cantilevered  truss  can  be  modelled  based  on 
the  Euler-Bernoulli  beam  formulation.  The  bending  motion  is  discribed  by  a  fourth 
order  partical  differential  equation  (PDE). 

pAQ  +  EIQ"''  =  0  (2.67) 


where  El  is  the  bending  stiffness  of  the  truss,  .4  is  the  cross-sectional  area  and  p  is 
the  density.  -The  independent  variable  along  the  length  of  the  truss  is  t  is  time, 
and  Q{zA)  is  the  flexural  displacement. 

The  boundary  conditions  (BC)  with  one  fixed  end  and  one  with  an  end  mass  are 
given  as 


g(o,a  =  o 

e'(o,o  =  0 


(2.68) 


To  solve  for  Q{z,t),  one  can  follow  the  standard  Fourier  method  and  convert  the  PDE 
into  infinite  assumed  modes  which  are  driven  parallelly  by  the  boundary  control  u{t). 
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Each  mode  can  be  expressed  as  a  second  order  underdamped  oscillation  equation. 


i^\ 


(2.69) 


(‘li  +  '2iiWiqi  +  w]qi  =  kiu{t),  i  =  1, 2,  •  •  • 


(2.70) 


where  (l)i{z)  is  the  spatially  dependent  mode  shape  coefficient,  qi{t)  is  the  time  depen¬ 
dent  modal  state,  tn,  is  the  modal  frequency,  if,-  is  the  damping  ratio,  and  ki  is  control 
gain  for  each  mode.  The  system  parameters  are  given  in  [19]  and  [20]. 

Control  action  is  obtained  using  proof-mass  actuators  with  neglected  actuator  dy¬ 
namics.  Sensing  of  truss  motion  is  accomplished  through  a  piezoelectric  accelerome¬ 
ter  mounted  on  the  structure  to  measure  the  relative  velocity  between  the  proof-mass 
actuator  and  the  structure. 

Combining  the  dynamic  equation  ((2.70))  with  the  actuation  and  sensing  mechanism, 
we  come  up  with  the  following  I/O  relationship: 

y{t)  =  yo{t)  +  '^!j,{t)  (2.71) 

yi{t)  +  2(iWiyi{t)  wfyi{t)  =  ku{t),i  =  1,2,  -  •  • ,  (2.72) 

where  y(t)  is  the  velocity  measurement  obtained  from  the  output  of  the  accelerometer; 
yo{t)  is  a  low  frequency  (~  0.1  Hz)  measurement  noise;  ;(/,(t),  i  =  1, 2,  •  •  •  is  the  flexural 
velocity  for  each  mode. 

With  the  system  parameters  Wi,  ki,  we  obtain  a.  state  space  representation  tor  each 
mode  of  interest 


Xi  =  AiXi  -t-  BiU, 

yi  =  CiXi,  (2.73) 

where  Xi  =  [q,  is  an  2x1  modal  state  vector,  u  is  th  e  scaler  control,  yi  is  the 
bending  velocity  of  mode  i, 


Ai  = 


0  1 

-wf  -2^iWi 


B,= 


bi 


,  Ci  -  [0  c,' 


(2.74) 


and 


biCi  —  ki. 


(2.75) 
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There  are  infinitely  many  possible  choices  for  bi  and  ci  satisfying  ((2.75)).  For  sim¬ 
plicity,  we  choose 

bi  =  \/m,  Ci  =  sgn{ki)^/\k\  (2.76) 

The  sampling  frequency  throughout  this  experiment  is  100  Hz,  which  corresponds  to 
sampling  period  T  =  0.01  second.  The  experimental  results  presented  here  are  for  the 
first  flextural  mode  only.  To  demonstrate  the  effectiveness  of  the  proposed  control 
method,  we  take  the  following  senario: 

Part  1;  0  <  t  <  2  Settlement  -  Uk  =  0. 

Part  2:  2  <  t  <  12  Shake  up  the  truss  -  Uk  =  0.3  sin  (toi  •  kT). 

Part  3:  12  <  f  <  20  Sliding  mode  control  v.s.  No  control. 


Figure  2.3:  Frequency  content:  Sliding  mode  control  (c=l)  v.s.  No  control 
The  sliding  surface  for  the  second  order  vibrating  system  is  chosen  as 
s(a;)  =  [c  Ijxi  =  CQi  +  qi  —  Q. 
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I 


Figure  2.4:  Frequency  content:  Sliding  mode  control  (c=10)  v.s.  No  control 

Experiments  are  conducted  for  c  =  1  and  c  10.  A  saturation  function  is  employed 
to  limit  the  applied  control 

Uk  =  Sat(ul\0d). 

Figure  2.3  depicts  the  frequency  content  (FFT)  of  y/,.  in  Part  3  for  c  =  1.  Similar 
result  for  c  =  10  is  shown  in  Figure  2.4,  in  vvliich  better  performance  is  obtained. 
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3.  VIBRATION  SUPPRESSION  IN  FLEXIBLE  STRUCTURES  VIA 
GENERALIZED  SLIDING  MODES 


In  this  chapter  the  problem  of  supressing  vibrations  in  flexible  sructures  by  means 
of  sliding-mode  control  is  considered.  Sliding  mode  control  algorithms  for  flexible 
structures  have  been  used  previously,  but  these  were  based  on  finite-dimensional 
models  (see  for  example  [22]  and  [23]). 

Sliding  mode  control  provides  stable  and  robust  closed  loop  systems  for  finite-dimensional 
plants,  but  the  direct  application  of  this  design  technique  based  on  approximate 
models  of  flexible  structures,  may  lead  to  undesirable  performance  of  the  infinite¬ 
dimensional  system.  This  phenomenon  is  due  to  the  delays  which  are  inherent  in 
such  systems.  The  presence  of  delays  can  be  explained  by  wave  propagation  across 
the  structure.  If  the  control  is  applied  to  the  boundary  of  the  flexible  structure,  its 
action  can  influence  other  parts  of  the  structure  only  after  the  wave  caused  by  the 
actuation  propagates  and  reaches  those  parts.  In  the  case  of  large  flexible  structures, 
these  delays  are  not  small  and  cannot  be  neglected.  The  application  of  sliding  mode 
control  ignoring  these  effects  leads  to  chattering  and  may  not  be  successful. 

Here  we  consider  mathematical  models  in  the  form  of  partial  differential  equations 
(PDE)  that  permit  us  to  take  into  account  the  features  described  above.  Therefore, 
more  appropriate  control  algorithms  may  be  designed. 

The  generalization  of  the  sliding  mode  control  concept  to  systems  with  delays  and  for 
more  general  dynamic  systems  described  by  semigroups  of  state  space  transformations 
was  originally  considered  in  [4],  [24].  We  base  our  approach  on  this  concept,  which 
can  be  called  manifold  control.  The  design  procedure  is  divided  into  two  steps.  In  the 
first  step,  the  manifold  is  designed  in  such  a  way  that  it  will  be  an  integral  manifold 
for  the  closed  loop  system.  In  the  second  step,  the  control  which  forces  the  system 
to  move  along  this  manifold  must  be  found.  The  crucial  point  of  the  problem  is  to 
design  the  integral  manifold  which  guarantees  system  stabilization. 

In  this  chapter,  the  nondispersive  wave  equation  is  chosen  as  a  canonical  form  for 
distributed  parameter  systems  described  by  partial  differential  equations.  Since  the 
nondispersive  wave  equation  is  equivalent  to  a  system  with  delay,  this  allows  for 
the  transformed  system  to  use  sliding  mode  control  algorithms  developed  earlier  for 
systems  of  differential-difference  equations  [4]. 

Subsequently,  a  class  of  flexible  structures  with  fourth  order  spatial  partial  derivatives 
is  considered  by  using  the  same  technique,  'hhis  is  accomplished  by  going  through 
a  two  stage  process  where  they  are  initially  Iraiislormed  into  the  second  order  form 
analyzed  originally. 
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The  problem  of  designing  the  control  law  which  assigns  the  desired  stable  integral 
manifold  to  the  system  can  be  solved  by  using  any  of  several  different  methods  (in¬ 
cluding  linear  methods).  The  use  of  sliding  modes,  however,  makes  the  closed  loop 
system  highly  insensitive  to  external  disturbances  and  parameter  variations. 

3.1  Generalization  of  the  Sliding  Mode  Control  Concept 

There  are  two  aspects  in  traditional  sliding  mode  control  design:  the  choice  of  the 
sliding  surface  and  control  synthesis  in  the  reduced  order  space.  From  the  point  of 
view  of  dynamic  system  theory,  the  sliding  surface  is  just  a  stable  integral  manifold 
of  the  closed  loop  .system,  with  a  specific  property  that  in  the  area  of  attraction  the 
system  state  is  absorbed  by  the  manifold  in  finite  time. 

For  finite  dimensional  systems  in  IR”  modeled  as 

(3.1) 

such  manifolds  can  exist  only  if  the  right  hand  side  does  not  satisfy  the  well  known 
Lipschitz  condition 

f{y)\<  L\x-tj\,  (3.2) 

which  is  usually  used  to  guarantee  the  uniqueness  of  the  solution  both  for  t  >  to  and 
t  <  to. 

Let  F{t;to,  Xo)  be  a  solution  of  the  system  of  ordinary  differential  equations  (3.1) 
with  initial  condition  x{to)  =  Xo,  i.e.  a  transition  function.  Then  the  Lipschitz 
condition  implies  that  F  is  defined  for  t  >  to  and  t  <  to-  The  family  of  state 
space  transformations  F{t',  to, .)  is  a  group  with  respect  to  the  composition  operation. 
The  inverse  element  to  F{t;to,.)  is  F{to;t,.).  For  an  asymptotically  stable  integral 
manifold,  the  trajectory  initiated  in  its  vicinity  tends  to,  but  never  reaches  it. 

In  contrast  to  systems  of  equations  whose  right-liand  sides  satisfy  the  Lipschitz  condi¬ 
tion,  in  systems  with  discontinuities  there  are  integral  manifolds  which  can  be  reached 
in  finite  time. 

Consider  the  system  in  IR" 

X  =  f{x)  +  B{x)u,  (3.3) 

where  f{x),  B{x)  are  functions  which  satisfy  tlie  Lipschitz  condition,  and  u  €  K™  is 
discontinuous  on  the  smooth  surfaces  {.r  ;  Si{x)  =  0}  i  =  1,2,. . .  ,m  in  IR": 
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(3.4) 


„  (r)  =  I  ^  ^ 

I  Uj.  (x)  if  5j(.t)  <  0 

If  the  sliding  mode  exists  on  the  intersection  of  discontinuity  surfaces  a  =  : 

Si{x)  =  0}  then 

1.  cr  is  an  integral  manifold  since  it  consists  of  the  system  trajectories; 

2.  the  uniqueness  of  the  inverse  to  the  shift  operator  F{t\to., .)  does  not  hold  on  a 
since  each  piont  on  cr  can  reached  from  at  least  two  different  points:  from  outside 
(because  of  the  finite  time  of  convergence)  and  from  points  on  the  manifold. 
This  is  a  characteristic  feature  of  sliding  mode  and  it  is  taken  as  a  basis  for  the 
generalization  of  the  sliding  mode  concept. 

The  general  definition  of  a  dynamic  system  in  any  metric  space  .T  (including  those 
under  consideration  described  by  partial  differential  equations)  utilizes  a  description 
in  the  form  of  a  transition  operator  F  similar  to  the  one  considered  above.  The 
following  Definition  was  introduced  in  [4] 


Definition  1  .r  6  X  is  said  to  be  a  sliding  point  at  the  time  instant  t  if  an  equation 
Fit',  ^0)0  =  0  for  every  to  <  t  has  more  than  one  .solution 

This  Definition  implies  that  the  sliding  manifolds  are  asymptotically  stable  manifolds 
to  which  the  system  state  converges  in  finite  time  from  any  initial  condition  in  the 
area  of  attraction. 

The  underlying  philosophy  of  the  proposed  approach  as  applied  to  PDE  models  of 
flexible  structures  is  the  same  as  in  any  sliding  mode  control  design  for  finite  di¬ 
mensional  systems.  After  representing  the  system  in  a  “convenient”  form,  a  sliding 
manifold  is  chosen  and  then  the  control  is  designed  such  that  the  system  state  reaches 
this  manifold  in  finite  time  and  then  “slides”  along  it.  The  control  in  this  case  is  not 
necessarily  discontinuous. 


3.2  Block  Control  of  Differential-Difference  Systems 

In  this  section  we  shall  consider  a  class  of  composite  systems  in  block  form  [25]  and 
[26]  comprized  of  blocks  of  differential  equations  coupled  with  blocks  of  difference 
equations.  We  shall  demonstrate  that  sliding  mode  controllers  may  be  effectively 
utilized  in  driving  the  state  to  zero. 
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We  consider  two  configurations: 
Configuration  A 


x{t)  —  Aiix{t)  +  Ai2z{t)  (3.5) 

z[t)  =  A2ix(t  —  t)  A22z{t  —  t)  BoU{t  —  t),  (3-6) 

Configuration  B 

z{t)  =  Auz{t  -  t)  +  Ai2x{t  -  t)  (3.7) 

x{t)  =  A2\z[t)  +  A22x{t)  +  Bou{t).  (3.8) 


It  is  assumed  that  x  G  IR"',  s  G  IR”^  and  u  G  IR’".  An,  A12,  A21,  A22,  Bo  are  matrices 
of  appropriate  dimensions.  Since  A12  is  not  necessarily  a  full  rank  matrix,  we  can  use 
the  representation: 


Ai2  =  5iC'2,  (3.9) 

where  Bi  and  C2  have  full  column  and  row  rank  respectively.  We  assume  that  the  pairs 
(All,  .^1)  in  both  types  of  configurations  are  controllable  and  the  system  ((72,  A22,  .^o) 
is  invertible. 

Denoting  v  =  C2Z  we  consider  it  as  a  new  control  variable  in  (3.5),  (3.7)  and  as  the 
output  variable  for  (3.6)  and  (3.8). 

In  the  first  block  of  Configuration  A 


x{t)  —  Aiia:(t)  +  B\v{t) 
we  use  sliding  mode  control  u(x)  =  co/(ui, . 


(3.10) 


,  e„2)  to  provide  stability 


Ui(.'r) 


_  (  v^{x)  if  Si(x)  >  0 
”  1  u~(a;)  -  if  Si{x)  <  0. 


(3.11) 


The  sliding  manifold  for  this  block  is 
ao  =  {x  :  s{x)  =  0}, 


(3.T2) 


where  s  =  col{si, . . . , Sn2)- 

Since  the  second  block  (3.6)  is  the  difference  system,  under  the  invertibility  condition 
it  is  possible  to  find  a  control  u  which  provides  the  desired  discontinuous  function  for 
the  output  v{x)  —  C2Z  [4].  This  equality  defines  another  manifold  in  the  system  state 
space 
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CTi  =  {(a;,r)  :  u(a;)  -  C2Z  =  0}. 


(3.13) 


As  a  result,  the  state  slides  on  the  intersection  a  =  cto fieri.  That  means  we  can  assign 
any  desired  rate  of  stability  for  the  closed  loop  system  (3.5)-(3.6). 

For  the  system  in  Configuration  B  the  design  procedure  is  different.  On  the  first 
step  we  find  a  linear  control  v  =  Dz{t)  for  the  first  block  (3.7) 

.::(<)  =  Anz{t  -  r)  +  Biv{t  -  r)  (3.14) 

to  establish  stability  for  the  difference  system  (3.7),  which  can  be  achieved  by  using 
the  sliding  mode  control  on  some  manifold  (72  system  as  described  in  the  previous 
section.  Then  the  problem  is  to  steer  s  =  D~{t)  —  C-ix-it)  to  zero  or  to  reach  the 
manifold 


(73  =  {(x,2:)  :  s(x,z)  =  0}.  (3.1.5) 

We  solve  this  problem  in  the  class  of  sliding-mode  control  algorithms: 

(3.16) 

The  resulting  motion  again  will  occur  on  the  intersection  cr  =  (72  f  (73. 


uf{x)  a  Si{x)  > 
u~{x)  if  5,(;r)  <  0 


3.3  Sliding-Mode  Control  for  Suppressing  Vibrations  of  a  Flexible  Bar 

As  a  first  example  we  consider  the  longitudinal  or  torsional  oscillations  of  a  flexible 
bar.  The  control  is  assumed  to  be  a  force  or  torque  applied  at  one  end  of  the  bar,  the 
other  end  is  free.  Let  Q  be  the  displacement  of  the  l)ar  from  the  unexcited  position. 
VVe  then  have  the  following  equations  for  a  unit  bar  with  normalized  parameters  [27]: 


d'^Q{t,x)  d^Q{t,x) 

(3.17) 

dB  dx"^ 

aQ(t,o) 

fir  = 

(3.18) 

o' 

II 

5 

(3.19) 
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where  x  is  the  position  along  the  rod  and  u{t)  denotes  the  actuation  force  or  torque. 
Applying  the  Laplace  transform  to  equation  (3.17)  and  boundary  conditions  (3.18), 
(3.19)  under  zero  initial  conditions 


g(0,.'c)  =  0,  (3.20) 

^Q(0,x-)  =  0,  (3.21) 

we  will  have: 

jy^Qip.x)  =^Q"ip,x)  (3.22) 

Q\p,0)  = -u{p)  (3.23) 

Q'{pA)  =  0,  (3.24) 


where  Q{p,x)  =  CQ{t,x),  u{p)  —  Cu{t). 

The  solution  of  this  boundary  value  problem  for  the  ordinary  differential  equation 
(3.22)  is 


Qip,^)  = - - 

gP  _  g-p  p 


(3.25) 


The  solution  of  the  stabilization  problem  depends  greatly  on  what  point  of  the  bar 
is  considered  as  the  system  output.  We  shall  consider  the  free  end  of  the  bar  as  an 
output  (thus  this  is  the  noncollocated  actuator/sensor  case). 

y{t)  =  Q{t,l).  ■  (3.26) 


From  (3.25)  we  obtain 


y{p)  =  <?(p,  1)  =  p  "  ^p-Hp)-  (3-27) 

eP  —  e  ^  p 

In  the  time  domain  the  correspondence  between  u{t)  and  y{t)  may  be  written  as 

y(t  +  l)  -  y{t  -  l)  =  2n{t)  (3.28) 
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or 


il{i)  -  y{t  -  2)  =  2u{t  -  1). (3.29) 

We  can  now  write  this  equation  in  the  form  of  the  differential-difference  system  Con¬ 
figuration  A  by  introducing  a  new  variable 

m  =  zit)  (3.30) 

z{t)  =  z{t-2)+2u{t-l).  (.3.31) 

The  block  representation  of  the  differential-difference  system  simplifies  the  develope- 
ment  of  the  control  algorithm.  Considering  the  variable  z{t)  in  the  first  block  (3.30) 
as  the  control,  we  can  obtain  sliding-mode  [1]  by  assigning 

s(t)  =  -Xsgn{rj{t)).  (3.32) 

This  ecjuality  is  valid  if 

s(t)  =  z(t)  +  Xsgn{y{t))  =  0.  (3.33) 

To  achieve  the  above,  we  can  use  the  control 

-  1)  -  \>^sgn[y{t  -|-  1)).  (3.34) 


This  control  algorithm  seems  to  be  noncausal,  however  using  an  extrapolator  it  can 
in  fact  be  realized  as  an  operator  on  the  current  and  the  past  values  of  the  con¬ 
trol  variables.  To  demonstrate  this,  solve  equation  (3.30)  taking  y{t)  as  an  initial 
condition; 

y{t  -1- 1)  =  y{t)  -1-  z{T)dT.  (3.35) 

or  using  (3.31) 

y{t+l)  =  y{t)+  [  {z{T-l)+2u{T))dT  =  y(t)+y{t-l)-y{t-2)+2  f  i((r)fir(3.36) 
Jt-1  Jt-1 

since  y{t)  =  z{t).  Substituting  y{t  +  1)  from  this  expression  to  (3.34)  and  again  using 
the  fact  that  z{t)  =  y{t)  we  will  have: 

u{t)  = -l-y{t  -  l) -l-Xsgn{y{t)  +  y{t  -  [)  -  y{t -2) +  2  (  if(r)dT).  (3.37) 

2  2  Jt-\ 
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With  this  control  the  system  (3. 30), (3. 31)  and  therefore  (3.17)  is  stabilized  in  finite 
time. 


3.4  Dispersive  Flexible  Structures:  an  Integral  Transformation  Approach 

The  problem  of  supressing  vibrations  in  the  systems  described  by  second  order  partial 
differential  ecjuations  as  studied  in  the  previous  sections  was  reduced  to  those  for  a 
differential-difference  system.  The  delay  in  these  systems  is  due  to  the  nondispersive 
properties  of  such  structures.  The  Euler-Bernoulli  beam  is  an  e.xample  of  a  dispersive 
structure.  In  this  section  we  consider  the  Euler-Bernoulli  beam  and  other  dispersive 
structures  of  the  second  and  fourth  order. 

Consider  now  the  problem  of  supressing  normal  vibrations  along  a  unit  length  flexible 
beam  described  by  equations  of  fourth  order.  One  end  of  the  beam  is  assumed  to  be 
clamped  while  a  control  force  is  applied  to  the  other.  The  Euler-Bernoulli  model  of 
the  beam  with  normalized  parameters  is  : 


d'^Qit^x)  d‘^Qit.x) 

(3.38) 

df^  'dx^ 

g(t,o)  =  o 

(3.39) 

Q'{t,o)  =  o 

(3.40) 

Q"(C1)  =  0 

(3.41) 

1)  =  u{t). 

(3.42) 

where  '  denotes  partial  differentiation  with  respect  to  the  spatial  variable  x. 

The  main  idea  of  our  approach  is  to  reduce  tlu'  order  of  the  controlled  part  of  the 
system  by  applying  an  integral  transformation: 

-P(C0=/  V{i,x)Q(i,x)dx.  (3.43) 

Jo 

Here  is  a  new  controlled  variable  and  .f  Is  a  new  independent  spatial  variable 

(0  <  ^  <  oo).  The  kernel  of  the  transformation  P  is  assumed  to  satisfy  the  same  type 
of  boundary  value  problem  as  Q  but  with  homogeneous  boundary  conditions. 
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_  crv{(,x) 
dx^  ■ 

(3.44) 

O 

II 

O 

(3.45) 

o 

II 

o 

(3.46) 

II 

o 

(3.47) 

v'"U,i)  =  0. 

(3.48) 

^  in  this  equation  is  analogous  to  a  time  variable  and  its  values  can  change  from  zero 
to  infinity. 

It  has  been  shown  [28]  that  under  these  conditions  F(t,<f)  satisfies  an  equation  of  the 
second  order  with  control  in  the  right  hand  side: 

cXP(M)  d‘P(U) 

gt2  =  gg;,  +'^(f)“(i), 

(3.49) 

where 

<^(0  = 1).  (3.50) 

In  this  equation  in  contrast  to  that  with  T>,  ^  is  a,  spatial  variable. 

We  can  say  that  the  transformation  (3.43)  “absorbs”  the  dispersive  properties  of 
the  fle.xible  structure  as  equation  (3.49)  describes  how  the  waves  are  travelling  there 
without  changing  their  form.  For  more  details  on  this  approach  and  a  more  general 
formulation,  see  [28]. 


3.5  Use  of  Neural  Networks  with  the  Integral  Transformation  Approach 

For  a  given  point,  ^  fhe  integral  transformation 

P{t,C)=  I' DiC,x)Q(t,  x)dx  (3.51) 

Jo 

may  be  realized  with  the  aid  of  an  artificial  neural  network  as  shown  in  Figure  3.1. 
As  shown  previously,  can  be  realized  as  the  output  of  a  differential-difference 
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system  of  the  form  of  (3.5)-(3.6)  with  input  u{t).  Rather  than  solving  for  D{^*,x)  and 
performing  the  integral  transformation  analytically,  a  neural  network  may  be  trained 
to  perform  these  tasks.  It  should  be  noted  that  the  transformation  is  a  functional 
depending  only  on  x  and  ^*-not  on  t. 


v(t) 


beam 


Neural  Network 


Model 


Figure  3.1:  Configuration  for  Training  Neural  Network 


The  input  u{t)  is  assumed  to  excite  both  the  beam  and  the  differential-difference 
model  (3.5)-(3.6).  The  output  of  the  neural  network  i/‘  =  0*  is  compared  to  = 
P{t,^*)  which  is  the  output  of  the  differential-difference  model.  The  error  (5‘  =  y^  —  y*' 
is  then  used  to  train  the  neural  network  using  some  training  algorithm  such  as  the 
standard  backpropagation  techniques.  [29,  30] 

The  activation  function  of  each  neuron  may  arbitrarily  chosen  to  be  any  differentiable 
function  At  some  time  t,  the  activation  of  the  neuron  in  any  of  the  layers 

is  given  by 

a;  =  X{h<]  (3.52) 

where  h\  is  the  total  input  for  the  neuron 

hi  =  Y^Wijd]  +  ei  (3.53) 

J 

where  dj  is  the  output  of  the  neuron  in  the  preceding  layer. 

Using  the  standard  quadratic  error  measure  for  time  t. 


2 


2 


(3.54) 


Vj-V 


t 


7 


and  the  overall  measure  of  error  E  =  "^E^  the  neural  network  is  trained  using  a 

t 

gradient  descent  algorithm  to  adjust  the  weights.  Variations  to  the  standard  back- 
propagation  algorithm  such  as  inclusion  of  a  momentum  term  in  the  update  rule  have 
been  found  to  be  useful  in  speeding  the  convergence  of  the  network 

dE 

Awij{t  +  1)  =  -7/7^ - h  aAwij{t) 

where  Wij{t  +  1)  =  Wij{t)  +  Awij{t). 

For  proper  training  of  the  neural  network,  all  modes  of  interest  of  the  beam  must 
be  excited.  Therefore,  the  input  to  the  beam,  u(t),  or  any  initial  conditions  must  be 
sufficiently  rich  to  excite  these  modes  of  interest.  Also,  due  to  the  computational  in¬ 
tensity  of  training  an  artificial  neural  network,  it  is  advantageous  to  train  the  network 
off-line. 

Once  the  neural  network  has  been  trained,  the  network  becomes  the  first  stage  of  the 
controller  as  shown  in  Figure  3.2.  The  output  of  the  neural  network  is  fed  to  the 
sliding- mode  controller  which  has  been  designed  for  the  differential-difference  model. 


(3.55) 


v(t) 


beam 


Q(t,x) 


Neural  Network 


Controller 


Figure  3.2:  Closed  Loop  System  with  Neural  Network 


Simulations  of  an  Euler- Bernoulli  beam  example  and  more  details  of  this  work  may 
be  found  in  [31]. 
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4.  RADIAL  BASIS  FUNCTION  NEURAL  NETWORKS  AND 
WAVELET  NEURAL  NETWORKS 


The  application  of  artificial  neural  networks  for  solving  learning  problems  and  for 
approximating  nonlinear  functions  has  occurred  in  various  fields  at  a  frenzied  pace. 
A  good  summary/introduction  into  the  field  of  neural  networks  can  l^e  found  in  [29]. 
When  solving  general  pragmatic  learning  problems,  the  training  data  tends  to  be 
irregularly  distributed  over  the  input  space.  The  radial  basis  function  methodology  is 
an  approach  for  multivariable  interpolation  which  is  especially  well-suited  for  dealing 
with  irregularly  positioned  data  points.  The  radial  basis  function  technique  constructs 
a  function  space  which  is  linear  and  depends  on  a  distance  measure  between  the  known 
data  points.  The  idea  of  a  multiresolution  learning  network  or  wavelet  network  is  to 
have  the  capability  of  dealing  with  irregularly  positioned  data  points,  generalizing  the 
nonlinear  relationship  between  the  given  inputs  and  desired  outputs,  and  being  able 
to  capture  the  desired  fine  details  of  the  nonlinear  mapping.  This  chapter  presents  a 
general  procedure  for  designing  neural  networks.  The  procedure  is  applicable  to  both 
radial  basis  function  neural  networks  and  wavelet  neural  networks. 


4.1  Radial  Basis  Function  Networks 

A  radial  basis  function  (RBF)  neural  network  (RBFNN)  is  a  network  composed  of 
input  vectors  X;  =  [  a,’,i  ;r,:2  •  •  •  Xip  ,  X,  G  IRP  and  output  vectors 

Y(  =  [  yn  yi2  •  ■  •  yiq  ]  ,  Y,-  G  for  i  =  1.2 _ ,r  exemplars  with  a  hidden  layer 

of  radially  symmetric  kernel  activation  functions.  .According  to  Broomhead  and  Lowe 
[32],  the  RBFNN  is  given  by  the  mapping  y,:/,  :  IW’  R,  Vi  =  1, 2, ... ,  /■ 

yik  =  Xok  +  E  1.2,...,y  (4.1) 

where:  Y;  =  [  y,i  y,'2  •  •  •  yiq  ] 

yik  is  Rh  output  element  for  ith  exemplar 

Aoi  represents  the  output  bias  for  the  Atli  element  of  the  output  vector 
L  is  the  number  of  kernel  neurons 

Xjk  is  the  weight  between  jth  kernel  node  and  kth.  output  element 
(f>{z)  :  =  Radial  Basis  Function 

Cj  G  represents  the  centroid  vector  for  tlie  jth  kernel  RBF 
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The  solution  to  (4.1)  yields  a  set  of  linear  equations  for  the  NN  weights  Xjk  which 
can  be  given  as  follows: 


Making  the  following  definitions: 

r  2/11  •••  .!/i? 

[t.  ■■  T,  ]  4  :  ..  (4.:!) 

i/ml  '  ’  ’  IJinq 

-^01  •  •  •  Xqij 

All  •  ■  •  Ai^ 

A-v/i  ■  •  •  Ail/, 

1  (/>(||Xi-C,||)  •••  ^(||Xi-Ca/||) 

1  <^(||X2-Ci||)  •••  ,^(||X2-Cm||) 

1  </>(||X„-Ci||)  •••  <^(||X^-Cv/||) 

and  closely  examining  (4.2)  shows  that  the  vectors  of  network  weights  Aj  can  be  solved 
for  independently  by  considering  only  the  Ti-  output.  Since  the  system  of  equations 
(4.2)  are  decoupled,  the  RBF  network  can  be  uniquely  determined  by  considering  it 
as  q  multiple  input,  single  output  (MISO)  systems.  The  NN  output  weights,  Xjk,  can 
be  found  by  solving  the  corresponding  linear  least  squares  optimization  problem. 

When  developing  a  RBFNN,  the  primary  design  parameters  to  be  determined  consist 
of  choosing  appropriate' centroids,  Cj,  and  smoothing  factor(s),  cTj.  RBFNN  perfor¬ 
mance  is  critically  dependent  upon  the  selected  centroids  and  smoothing  factors.  Sev¬ 
eral  methods  are  commonly  used  for  selecting  appropriate  centroids.  One  approach 
uses  the  standard  k-means  clustering  algorithm  [33],  [34],  as  an  iterative  process  over 
the  entire  set  of  training  data.  The  other  approach  uses  the  standard  k-means  clus¬ 
tering  algorithm  as  an  adaptive  process  in  real  time  [34].  The  common  method  for 
determining  the  smoothing  factor(s)  is  the  heuristic  known  as  “P  nearest- neighbor’' 
[35].  The  primary  problem  encountered  wuth  these  selection  procedures  is  that  the 
resultant  RBF  networks  tend  to  have  a  large  number  of  kernel  nodes,  and  some  even 
have  poor  performance.  Chen,  Cowan  and  Grant  [36]  have  developed  a  systematic 
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approach  to  selecting  the  most  significant  centers  from  a  set,  usually  large,  of  can¬ 
didate  centers.  The  approach  by  Chen,  et  al.  [36]  uses  an  orthogonal  least  squares 
(OLS)  procedure  in  a  forward  regression  manner  to  determine  the  significant  centers 
to  within  a  specified  tolerance  p,  0  <  p  <  1.  The  selection  of  the  tolerance  p  is  the 
key  to  balancing  the  accuracy  and  complexity  of  the  final  RBF  neural  network  design. 


4.2  Wavelet  Neural  Network  Architecture  1 

One  wavelet  network  structure  follows  readily  from  a  reconstruction  formula  for  any 
function  in  Hilbert  space,  f  E  H,  based  on  the  discrete  wavelet  transform  coefficients 
for  the  function  given  by  Daubechies  [37] 

/  =  A  1  D  i>m,n  /or(B/A-i)«l  (4.6) 

as  long  as  there  exist  A  >  0  and  B  <  oo  satisfying 

^ll/lt<  E  i{/.<A..,->t<-S||/ir  (4.7) 

m,n£2Z 


where:  {x,y)  represents  the  inner  product 

(/,  ■i/’m.n)  represents  the  discrete  wavelet  transform  coefficients 

m,n  {x)  =  ao'  —  nbo)  V???,  7Z  ^  2Z}  represents  a  discrete  family  of  wavelets 

A  and  B  represent  the  frame  bounds 


When  the  frame  bounds  satisfy  the  condition,  A  =  the  frame  becomes  a  tight 
frame,  and  if  =  1  \Im^n  E  4^,  then  the  ibra.n  constitute  an  orthonormal  basis 

for  7Y,  Daubechies  [37].  In  general,  frames,  whether  tight  or  not,  are  not  orthonormal 
bases,  but  instead  simply  span  7^,  i.e.  admit  redundancy  in  representing  any  /  G  77. 
Thus,  the  wavelet  network  given  by  equation  (4.8)  can  approximate  any  function  in 
77,  or  for  more  practical  purposes  C^[1R). 

One  form  of  the  feedforward  wavelet  neural  network  (WNN)  architecture  is  a  network 
composed  of  input  vectors  =  [  xn  Xi2  •  •  •  x,.p  ]^  ,  Xj  G  IR^  and  output  vectors 
Yi  =  [  yn  yi2  •  ’  •  yiq  ]  ,  Y*  G  for  i  =  U2,. . .  ,r  exemplars  with  a  hidden  layer 
of  wavelet  activation  functions.  According  to  [38], [39],  the  feedforward  WNN  can  be 
represented  by  the  mapping  yik  \  IRF  IR,  V/'  =  1/2, . . . ,  r 

L 

=  Aoi 0(f/^a5f(aj)X.  -  bj)  VA- =  1, 2, . . . ,  9  (4.8) 


where:  ^/’(z)  =  Mother  Wavelet  Function 

diag(aj)  represents  a  diagonal  matrix  whose  elements  are 

aj  =  diag{aO~"^J'^  G  BV’  is  the  dilation  vector  for  the  ;th  wavelet 

111  =  . . . ,  is  dilation  integer  vector 

bj  =  6on  is  the  translation  vector  for  the  ,yth  wavelet 

11  =  [uj,!, . . . ,  Uj^p]'^  is  translation  integer  vector 

Since  the  WNN  equation  (4.8)  has  the  identical  format  as  (4.1),  it  follows  that  the 
WNN  can  be  uniquely  determined  by  considering  it  as  q  MISO  systems. 


4.3  Multiresolution  Signal  Analysis 

The  theory  of  multiresolution  signal  decomposition  of  Mallat  [40]  provides  the  crucial 
framework  for  formulating  and  designing  neural  networks  with  wavelet  activation 
functions.  A  recapitulation  of  these  concepts  as  given  in  [40],  [41],  and  [37]  will  now 
be  presented. 

A  multiresolution  signal  analysis  entails  a  sequence  of  successive  closed  approximation 
subspaces,  Vm  with  m  G  2Z,  satisfying  the  following  properties: 


1. 

•  •  •  C  V^2  C  El  C  Eo  C  VTi  C  VT2  C  •  •  ■ 

(4.9) 

2. 

[j  Vm  =  C\M) 

(4.10) 

3. 

II 

(4.11) 

4. 

fix)  eVm<=>  fi2x)  G  En-1 

(4.12) 

5. 

f(x)GVo=^fix-j)eVo  Vjez 

(4.13) 

which  due  to  (4.12)  implies  f{x)  G  =>  j\x  —  2'”j)  G  Vm  Vj  G  Z 


6.  There  exists  a  scaling  function  </>  G  V'o  such  that 

j  £  is  an  orthonormal  basis  in  Vo  (4'14) 

with  =  2"’"/^^^(2"’”a;  -  j)  \/j,  )n  G  ^ 

Equations  (4.12)  and  (4.14)  imply  that  {4>,n.j\  j  G  is  an  orthonormal  basis  for 
Vm  Vm  G  Z.  Define  to  be  the  orthogonal  projection  operator  onto  Vm-  By  (4.10), 
limTO=>.-cc.  Pmf  =  f  V/  G  The  crux  of  multiresolution  analysis  is  that  if  a  set 
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of  closed  subspaces  satisfies  (4.9)  -  (4.14),  then  there  exists  an  orthonormal  wavelet 
basis  of  C^{M)  given  by: 


>^1,71  e  Z) 

t/Vn(x)-2-™/2i/^(2-"^.r-n)  (4.15) 

which  for  all  /  €  C^{M)  satisfies  the  relation 

Pm-lf  =  Pmf  +  Y,  {f^^m,n)i'm.n  (4.16) 

For  every  m  €  Z,  define  the  orthogonal  complement  of  Vm  in  Vm-\  as  Wm-  At  every 
scale 


Vm.l  =  VmeWm  (4.17) 

Wm  -L  Wm'  if  mj^m'  (4.18) 

Wm  C  Vm>  if  m  >  m'  (4.19) 

Then,  for  m  <  M, 

M  —  m  —  1 

Vm  =  VM(B  0  Wm-^  (4.20) 

t=0 

From  equations  (4.10)  and  (4.11),  it  follows  that  C^{M)  can  be  decomposed  into 
mutually  orthogonal  subspaces  according  to 

^Wm=^C^R)  (4.21) 


Lastly,  the  Wm  subspaces  acquire  the  scaling  property  of  equation  (4.12): 

fix)  G  Wm  f(2x)  G  (4.22) 

With  the  above  definition  for  Wm,  equation  (1.  Hi)  (‘(luivalently  states  that  {t/’m.n;  n7,  »  G 
Z]  is  an  orthonormal  basis  for  Wm  for  fixed  ni.  In  addition,  equation  (4.22)  guar¬ 
antees  that  if  {Tpoy,  n  G  Z]  is  an  orthonormal  l)a.sis  for  Ifb,  then  {f>m,n\  7n,7i  Z} 
will  also  be  an  orthonormal  basis  for  Wm,  for  any  in. 

The  orthonormal  basis  constructed  above  using  tlie  multiresolution  signal  analysis, 
given  by  equation  (4.15),  is  only  valid  for  C^{IR).  I'lie  easiest  way  of  constructing  an 
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orthonormal  basis  for  C^{W’)  is  to  take  the  tensor  product  functions  generated  by  p 
one  dimensional  wavelet  bases; 


mi,ni  ;m2,n2;'-;mp,np\^l  t  ^2i  ^  t  ^p)  — 
0  m  1  ,ni  (.•ri)'0 

m2  ,n2  (•^2)  •  •  •  i’mp,np{-rp) 


(4.23) 


where  {'^mi,ny,m2,n2\--,mp,np',  rui^Ui  G  Z,]  is  an  orthonormal  basis  for  C'^{1R!’).  An¬ 
other  construction  for  multidimensional  w'avelets  comes  from  performing  a  multidi¬ 
mensional  multiresolution  analysis  in  which  the  dilations  of  the  resultant  multidimen¬ 
sional  wavelets  control  the  p  variables  simultaneously  [37].  However,  multidimensional 
wavelets  generated  using  either  construction  suffer  from  the  drawback  of  depending 
on  the  dimension  p  of  the  input  space  [37].  The  generation  of  nonseparable  multidi¬ 
mensional  wavelets  appears  to  be  a  viable  technique  for  overcoming  the  dependence 
of  multidimensional  wavelets  on  the  dimension  p  of  the  input  space  [42]. 


4.4  Wavelet  Neural  Network  Architecture  2 

A  second  architecture  for  the  wavelet  network  comes  directly  from  a  multiresolution 
signal  analysis.  Let  /,„_i(a;)  be  the  appro.ximatiou  of  f{x)  at  the  (m  —  1)  scale.  Then, 
since  ^  Kn-i  and  recalling  equation  (4.17),  it  follows  from  equation  (4.16) 

that, 


./"m  — 1(^)  —  T  ^  ^  V’m.n (a' )  (4-24) 

If  f{x)  ^  /„^_l(x),  i.e.  the  (  m  —  1)  scale  corresponds  to  the  sampling  rate  of  the 
function  to  be  approximated,  and  recalling  ecpiation  (4.14),  then  equation  (4.24) 
becomes 

M 

y(ir)  ~  ^  ^  CMn0A/,n(^)  ^  ^  ^  ^  )  (4.25) 

-  i=m 

Note  from  equations  (4.14)  and  (4.12)  that  f.\i(x)  =  Y.aez  The  equation 

(4.25)  was  used  in  [43]  to  represent  a  wavelet  neural  network. 

Now  it  should  be  clear  that  the  wavelet  neural  network  given  by  equation  (4.8)  can 
alternatively  be  represented  by  the  mapping  /yn-  :  Ui'^  — >•  iR,  Vf  =  1,2,...,  r 

M  ;V 

ya-  =  Aoi+l]A,i  <i>{diag{a,)X,-h,)+Y^  Xp.  c(r/m^/(a,  )X.-b,)  V/:  =  1, 2, . . . ,  </(4.2()) 

/=!  J=l 
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4.5  Wavelet  Neural  Networks 


When  developing  a  WNN,  the  primary  design  parameters  to  be  determined  consist 
of  choosing  appropriate  translation  vectors,  6on,  and  dilation  factors,  This  is 

equivalent  to  selecting  kernel  wavelets  from  the  infinite  set  of  candidate  wavelets 
represented  by 

{0m,n(z)  =  {detAy^'^ip{Az  -  6on)m,n  €  e  M’  (4.27) 

where:  A  =  , . . . ,  aO~'"p).  For  multidimensional  input  spaces,  there  is 

an  additional  consideration  of  deciding  how  the  multidimensional  wavelets  are  to 
be  generated.  The  separable  multidimensional  wavelets  suffer  from  the  drawback  of 
being  dependent  on  the  dimension  p  of  the  input  space.  Another  approach  is  to 
generate  nonseparable  multidimensional  wavelets  [42].  Alternately,  a  radial  wavelet 
function  can  be  used,  which  is  defined  by  the  mapping  T  :  — )■  M  [39]: 

y{z)  =  ^(||z||2)  (4.28) 

Radial  wavelets  are  not  dependent  on  the  dimension  of  the  input  space  and  are 
therefore  better  suited  for  large  dimensional  learning  problems  than  direct  product 
wavelets.  Note  that  with  radial  wavelets  the  same  dilation  factor  is  used  over  each 
dimension  of  the  input  space,  i.e.,  diag[eij)  =  diag{aO~^^  ^ . . . ,  ). 

It  has  been  shown  that  for  appropriately  chosen  0,  uq,  and  bo,  equation  (4.27)  comprise 
frames  of  C^(IR^)  [39].  The  first  step  in  selecting  appropriate  kernel  wavelets  is  to 
truncate  equation  (4.27)  into  a  finite  set,  which  requires  a  pHoH  knowledge  about  the 
given  approximation  problem.  When  0  is  chosen  to  be  a  radial  wavelet,  the  truncated 
wavelet  set  can  be  denoted  by 

{^m,,n,(z)  =  ’"'^'^(ar'z  -  boiy)  m,  G  Z  C  n,  G  Z  C  (4.29) 

Any  truncated  wavelet  set  should  sufficiently  cover  the  domain  of  the  input  vectors 
at  the  desired  resolution  level,  rrij,  while  maintaining  the  position  of  the  wavelets, 
iij,  somewhere  within  the  domain  of  the  input  vectors  and  allowing  the  support 
of  the  wavelets  to  extend  beyond  the  range  of  the  input  training  data.  The  trun¬ 
cated  wavelet  set,  equation  (4.29),  has  the  form  of  a  regular  pyramid,  see  Figure  4.2, 
meaning  that  wavelets  corresponding  to  the  same  dilation  parameter  are  equally 
spaced/ translated  within  the  domain,  and  as  t  he  dilation  parameter  increases,  the 
corresponding  wavelets  become  denser/narrower  within  the  domain.  The  parallels 
between  the  RBFNN  and  the  WNN  should  now  be  apparent.  Choosing  the  centroids, 
Cj,  for  the  RBFNN  is  equivalent  to  selecting  the  translation  parameter  vectors,  ^onj. 
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for  the  WNN.  Similarly,  choosing  the  smoothing  factor(s),  cr^,  for  the  RBFNN  is 
equivalent  to  selecting  the  dilation  parameter(s),  ,  for  the  WNN.  Thus,  the  OLS 
procedure  [36]  can  be  used  on  the  truncated  set  of  candidate  wavelets  equation  (4.29) 
to  determine  the  significant  wavelets  to  within  a  specified  tolerance  p,  0  <  p  <  1, 
[44].  The  selection  of  the  tolerance  p  is  also  the  key  to  balancing  the  accuracy  and 
complexity  of  the  final  WNN  design. 


4.6  Network  Design 

Figure  4.1  below  shows  the  general  architecture  for  the  RBFNN  and  for  the  WNN. 


yi 


Xoi  ■  ^02 


^J1 


y: 

Output 
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.r''  ■ 
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Figure  4.1:  RBF  and  Wavelet  NN  Architecture  for  Learning 


Once  the  training  data  has  been  determined,  then  the  NN  designs  can  be  carried  out. 
The  kernel  node  activation  function  for  the  RBFNN  (4.1)  and  the  mother  scaling 
activation  function  for  the  WNN  when  using  the  second  architecture  was  selected  to 
be  a  multivariate  Gaussian  function. 


(f>{z)  =  exp  z  E  IR’’ 


(4.30) 


The  mother  wavelet  activation  function  for  both  architectures  of  the  WNN  (4.8) 
were  selected  to  be  the  Laplacian  of  the  multivariate  Gaussian  function,  i.e.  a  radial 
wavelet. 


V’(z) 


^4z^z  —  2p 


exp 


;  ze  R’’ 


(4.31) 
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where  <;  represents  the  variance  of  the  multivariate  Gaussian.  The  development  of 
the  NNs  involved  a  tradeoff  between  the  complexity  of  the  NNs  and  the  resultant 
performance  or  accuracy  of  the  NN  interpolation. 

Design  Goal:  Minimize  the  complexity  of  the  NNs  by  keeping  the  number  of  kernel 
nodes  to  a  minimum,  while  maintaining  the  normalized  error  below  5%. 

For  the  RBFNN  design,  since  a  RBFNN  with  the  same  smoothing  factor  in  each 
kernel  node  is  capable  of  universal  approximation  [45],  this  characteristic  should  be 
exploited  in  order  to  simplify  the  number  of  design  parameters  to  be  determined.  The 
initial  smoothing  factor  can  then  be  determined  using  a  nearest  neighbor  criterion. 
Now  the  critical  design  parameters  to  be  determined  are  the  locations  of  the  centroids 
for  the  kernel  functions. 

For  the  WNN  design,  a  grid  of  the  wavelets  at  the  multiple  I’esolution  levels  for  each 
input  range  must  be  determined  as  shown  in  Figure  4.2.  The  wavelet  coefficients  ciq 
and  bo  affect  the  alignment  of  the  grid  of  wavelets  for  each  input  range.  Since  for 
appropriately  chosen  ao  >  1  and  bo  >  0,  radial  wavelets  comprise  frames  of  C^{1RF) 
[39],  it  was  decided  to  utilize  this  property  in  order  to  simplify  the  number  of  design 
parameters  to  be  determined. 

Resolution 
Level ,  m 

f 

m  f  X 

coarse  \ 

•  » 

•  I  - 

"^fine  r  xxxxxxxxxxx 

j  [  Input  Range  1 

Wavelet  Position  ,  n 

Figure  4.2:  Grid  of  wavelets  at  multiple  resolutions 

In  order  to  form  the  grid  of  wavelets  for  each  input  range,  the  minimum  distance 
between  data  points  for  each  input  range  needs  to  be  determined,  Then  the 

finest  resolution  necessary  for  the  grid  over  each  input  range  is  given  by  ao  — 
which  yields: 

rriijine  ^  rOUnd{\og^Q  ^Xi^rnin)  (4.32) 

The  coarsest  resolution  for  the  grid  over  each  input  range,  nii^coarse:  is  determined  by 
the  largest  value  of  m  for  which  the  first  occurence  of  the  smallest,  nonzero  number 


of  wavelets  in  the  input  range  exists.  Now  the  set  of  dilation  integers  for  each  input 
range  is  given  by 


{^T'i,finei^i,fine  d”  f i  i  ^^i,ccarse}  —  C  ^  (4.t33) 

Once  the  corresponding  grid  for  each  input  range  has  been  determined,  then  the 
multidimensional  wavelets  need  to  be  generated.  The  use  of  radial  wavelets  greatly 
simplifies  this  task,  since  only  those  values  of  m  common  to  all  input  ranges  need  to 
be  considered.  Recall  that  the  same  dilation  factor  is  used  over  each  dimension  of  the 
input  space,  diag{sLj)  =  diag(aO~^"'i , . .  • ,  aO“'"J ),  therefore 

_  V  _ 

Z  =  Pi  Zj  ;  with  m-j  ^  Z  d  Z  (4.34) 

%■=! 

Recall  that  radial  wavelets  have  the  added  benefit  of  being  independent  of  the  dimen¬ 
sion  of  the  input  space.  Now,  the  corresponding  set  of  translation  integers  for  each 
common  resolution  level,  mj,  is  determined  by 

Nm,  =  {round{Xi/{aO^^JbO));i  =  I . r} 

Z  =  {Nnx^  '\/mj  eZcZ} 

ly  e  Nm,  c  Z  c  (4.35) 

Note  that,  for  each  m^,  there  are  |  Nm^  |,  (number  of  elements  in  the  set  Nmp, 
corresponding  Uj  vectors.  Then  the  truncated  wavelet  set,  equation  (4.29),  is  formed 
by  using  equations  (4.34)  and  (4.35). 

The  design  procedure  utilized  in  developing  RBF  and  wavelet  NNs  can  be  summarized 
for  a  given  set  of  training  data  as  follows: 

la.  Select  RBF  function,  ^(z).  OR 

lb.  Select  mother  wavelet,  t/’(z),  the  scaling  function  c)(z)  if  using  the  .second 
WNN  architecture,  and  the  dilation  constant,  oq. 

2a.  Select  the  RBF  smoothing  factor,  a.  OR 


2b.  Select  the  wavelet  ti'anslation  constant,  l\). 
3a.  Initialize  the  RBF  centroids,  Cj  =  X,.  OR 
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3b.  Determine  the  set  of  wavelet  dilation  integers,  G  Z,  by  generating  a 
grid  for  each  input  range  and  then,  for  each  iiij,  find  the  set  of  translation 
integers,  iij  G  Nm^  C  Z. 

4.  Choose  the  selection  tolerance,  p  (typical  0.0001),  and  use  OLS  algorithm 
[36]  to  determine  significant  RBF  centroids  or  wavelets  from  the  truncated 
set. 

5.  Get  least  squares  solution  for  weights,  Xjk. 

6.  Repeat  above  steps  as  necessary. 

sectionNeural  Network  Modeling  for  Control  Applications  A  RBFNN  and  a  VVNN 
will  both  be  developed  in  order  to  model  several  actuator  nonlinearities  in  a  flexible 
structure.  The  purpose  of  the  NNs  will  be  to  predict  the  dynamical  characteristics 
of  the  actuator  nonlinearity.  A  considerable  amount  of  uncertainty  exists  in  typical 
actuator  noniinearites.  The  output  vector,  Yi  —  [  ijn  yi2  ],  for  the  NNs  consisted 
of: 


1.  Desired  actuator  torque. 

Each  training  run  was  sampled  at  a  rate  of  lOHz  for  500  seconds,  resulting  in  5001  data 
points.  It  was  determined  that  a  minimum  of  6  training  runs  were  necessary  to  evenly 
partition  the  input  space,  and  these  6  training  runs  were  selected  as  the  preliminary 
set  of  training  data.  This  preliminary  set  of  training  data  contained  30,006  data 
points.  Since  the  kernel  node  centroids  are  initially  selected  to  correspond  with  the 
set  of  training  data,  the  centroid  selection  procedure  would  have  to  be  carried  out  over 
an  initial  number  of  30,006  kernel  nodes.  As  can  be  observed,  the  amount  of  computer 
memory  and  computational  time  necessary  for  even  the  minimal  set  of  training  data 
is  exorbitant!  Also,  if  additional  training  runs  are  deemed  necessary,  the  size  of  the 
set  of  training  data  wilf  only  increase.  One  solution  to  the  excessive  data  problem  is 
to  sample  the  training  runs  at  a  slower  rate  to  reduce  the  number  of  data  points  down 
to  a  manageable  number,  and  then  test  the  performance  of  the  resultant  NN  with 
the  lOHz  set  of  training  data.  In  theory,  as  long  as  the  training  data  is  chosen  to  be 
persistently  exciting  over  the  entire  input  space,  retraining  of  the  NNs,  i.e.  changing 
the  kernel  node  to  network  output  weights,  will  be  unnecessary.  In  other  words,  the 
resultant  NN  performance  with  non-training  set  data  should  be  comparable  with  its 
performance  over  the  set  of  training  data.  Additional  testing  data  not  included  in 
the  training  set  was  then  used  to  measure  the  performance  of  the  NNs. 
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4.7  RBF  Neural  Network  Simulation  Results 


The  RBFNN  designs  presented  in  this  section  follow  the  procedure  outlined  in  Section 
4.6  with  only  the  smoothing  factor,  a,  being  modified  between  designs.  The  initial 
smoothing  factor  for  each  kernel  node  w'as  determined  from  the  plot  of  the  nearest 
neighbor  distance  versus  training  data  vector  index.  The  mean  nearest  neighbor 
distance  is  0.98.  Based  on  the  nearest  neighbor  data,  the  initial  smoothing  factor 
was  selected  to  be  1.0.  Then  the  centroid  selection  procedure  [36]  was  employed  for 
a  tolerance  p  of  0.0001.  The  selection  procedure  resulted  in  260  kernel  nodes  being 
selected  and  an  output  bias.  The  main  trends  observed  in  the  RBF  designs  were  that 
the  accuracy  of  the  RBFNN  near  points  of  discontinuity  decreases  as  a  increases  and 
that  the  complexity  of  the  resultant  RBFNN  tends  to  decrease  as  a  increases. 


Scatter  Plot  of  Nearest  Neighbor  Distance  for  Neural  Network  Training  Data 


Figure  4.3:  Scatter  Plot  of  Nearest  Neighbor  Distance  for  NN  Training  Data 

Table  4.1  summarizes  the  RBFNN  designs  discussed  above  in  terms  of  the  design 
parameter  cr,  the  number  of  kernel  nodes  L,  and  the  normalized  error  for  each  training 
vector  of  exemplars,  Ek\  given  by  equation  ((4.36))  for  the  set  of  training  vectors. 


i=\ 


(4.36) 


where:  represents  the  normalized  error 

for  the  kth  training  vector 
r  is  the  number  of  exemplars 
in  each  training  vector 
t;  is  desired  output  for  fth  exemplar 
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Yi  is  the  NN  output  for  ith  exemplar 


a 

L 

El 

E2 

Es 

^4 

E, 

Ee 

1.0 

260 

0.81 

0.70 

0.40 

0.58 

0.26 

0.32 

5.0 

61 

0.45 

0.52 

0.59 

0.67 

0.70 

0.12 

7.5 

45 

0.47 

0.59 

0.55 

0.55 

0.58 

0.25 

10 

44 

0.70 

0.62 

0.57 

0.61 

0.38 

0.16 

15 

53 

0.71 

0.73 

0.51 

0.58 

0.32 

0.13 

25 

41 

0.82 

0.99 

0.82 

0.78 

0.39 

0.06 

Table  4.1:  RBFNN  Training  Design  Summary 


Table  4.2  compares  designs  3  and  4  when  tested  using  additional  data  not  included 
in  the  training  set. 


a 

L 

El 

E2 

Es 

E4 

Es 

Ee 

7.5 

45 

0.48 

1.17 

0.55 

1.07 

0.58 

0.25 

10 

44 

0.71 

1.24 

0.57 

1.11 

0.38 

0.16 

Table  4.2;  RBFNN  Testing  Design  Summary 


Remark:  If  the  normalized  error  was  5%  at  every  point  in  a  training  run,  =  250. 
Therefore,  the  results  shown  in  Tables  4.1  and  4.2  are  excellent  and  more  than  satisfy 
the  design  goal  stated  in  Section  4.6 

The  RBFNN  output  responses  are  given  in  Figures  4. 4-4. 6.  Figures  4. 4-4. 6  compare 
the  RBF  output  responses  for  design  4  with  the  actual  set  of  testing  data  and  provide 
the  resultant  error  and  normalized  error  responses.  Observation  of  Figures  4. 4-4. 6 
show  that  the  performance  is  excellent  and  the  goal  of  maintaining  the  normalized 
error  below  5%  everywhere  except  at  an  occasional  point  of  discontinuity  is  satisfied. 


4.8  Wavelet  Neural  Network  Simulation  Results 

The  WNN  designs  follow  the  procedure  outlined  in  the  Network  Design  Section  with 
the  translation  constant,  6o,  and  the  mother  wavelet  activation  function,  ec[uation 
((4.31)),  being  modified  between  designs.  The  mother  wavelet  was  modified  by  vary¬ 
ing  the  variance  parameter,  s.  Initially,  c  =  1.  The  dilation  constant  was  set  to 
ao  =  2  and  the  initial  translation  constant  was  chosen  as  6o  =  1.  The  resulting  grid  of 
wavelets  for  each  input  range  is  summarized  below  in  Table  4.3,  where  |  N  \  denotes 
the  number  of  elements  in  the  set  Nm  . 
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X4 

mi 

|iV| 

m2 

\N\ 

m3 

\N\ 

m4 

|iV| 

n 
-  ( 

1450 

-7 

1632 

0 

156 

-7 

52 

-6 

731 

-6 

832 

1 

78 

-6 

27 

-5 

366 

-5 

411 

2 

39 

-5 

14 

-4 

183 

-4 

207 

3 

21 

-4 

7 

-3 

92 

-3 

103 

4 

11 

-3 

4 

-2 

46 

-2 

53 

5 

5 

-2 

3 

-1 

24 

-1 

27 

6 

4 

-1 

2 

0 

12 

0 

13 

7 

2 

0 

1 

1 

7 

1 

7 

8 

1 

2 

4 

2 

5 

3 

3 

3 

3 

4 

1 

4 

1 

Table  4.3:  Grid  of  Wavelets  for  Actuator  Data 


From  Table  4.3  it  follows  that  Z  =  {0}  and  therefore  —  0.  The  set  of  translation 
integers,  G  were  determined  using  equation  ((4.35)).  Then  the  selection 

procedure  [36]  was  employed  for  a  tolerance  p  of  0.0001.  The  selection  procedure 
resulted  in  88  wavelets  being  selected  and  an  output  bias.  The  main  trends  observed 
in  the  wavelet  designs  were  that  the  accuracy  ot  the  WNN  increases  as  6o  decreases 
and  that  the  complexity  of  the  resultant  WNN  tends  to  decrease  as  s  increases. 

Table  4.4  summarizes  the  WNN  designs  discussed  above  in  terms  of  the  design  pa¬ 
rameters  6o  and  the  number  of  kernel  nodes  L,  and  the  normalized  error  for  each 
training  vector  of  exemplars,  Ek-,  given  by  eciuation  ((4.36)). 


bo 

L 

El 

E2 

Eo 

El 

E, 

Ee 

1.0 

1 

88 

11.4 

5.72 

0S.2 

50.1 

13.5 

1.35 

0.5 

1 

245 

0.81 

0.58 

0.53 

0.70 

0.34 

0.32 

0.1 

1 

242 

0.71 

0.62 

0.51 

0.76 

0.48 

0.39 

0.5 

4 

104 

0.46 

0.63 

0.45 

0.65 

0.94 

0.10 

0.5 

10- 

48 

0.51 

0.74 

0.41 

0.62 

0.52 

0.18 

0.1 

10 

46 

0.55 

0.61 

0.45 

0.66 

0.54 

0.24 

Table  4.4:  WNN  Training  Design  Summary 


Table  4.5  compares  designs  5  and  6  when  testc'd  using  additional  data  not  included 
in  the  training  set. 

The  results  shown  in  Tables  4.4  and  4.5  are  alsc;  (excellent  and  more  than  satisfy  the 
design  goal.  Comparison  of  Tables  4.2  and  4.5  demonstrates  that  the  performance  ot 
the  RBFNN  and  WNN  are  basically  identical  for  tliis  application.  The  WNN  output 
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bo 

L 

Ex 

E2 

Ez 

E4 

Es 

Ee 

0.5 

10 

48 

0.52 

1.35 

0.40 

1.05 

0.52 

0.18 

0.1 

10 

46 

0.56 

1.20 

0.45 

1.08 

0.54 

0.24 

Table  4.5:  WNN  Testing  Design  Summary 


responses  are  given  in  Figures  4. 7-4.9.  Figures  4. 7-4. 9  compare  the  WNN  output 
responses  for  design  6  with  the  actual  set  of  testing  data  and  provides  the  resultant 
error  and  normalized  error  responses.  Observation  of  Figures  4. 7-4. 9  shows  that  the 
performance  is  excellent  and  the  goal  of  maintaining  the  normalized  error  below  5% 
everywhere  except  at  an  occasional  point  of  discontinuity  is  satisfied. 


Neural  Network  Response  vs.  Actual  Response  for  Sigma  =  10  and  Rho  =  0.0001 


Erfor  in  Neural  Network  Response 


Normalized  Error 


time  (seconds)  time  (seconds) 


Figure  4.4:  RBFNN  Design  4  -  Testing  Run  1 
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Neural  Network  Response  vs.  Actual  Response  for  Sigma  =  10  and  Rho  =:  0.0001 


time  (seconds) 


time  (seconds) 


Figure  4.5:  RBFNN  Design  4  -  Testing  Run  2 


Neural  Network  Response  vs.  Actual  Response  for  Sigma  =  10  and  Rho  =  0.0001 


time  (seconds)  time  (seconds) 


Figure  4.6:  RBFNN  Design  4  -  Testing  Run  4 
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bc1{t)  ^  e(t)  bc1{t) 


Neural  Network  Response  vs.  Actual  Response  for  a0=2,b0=0.1  ,Rho=0.0001 
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Figure  4.7:  WNN  Design  6  -  Testing  Run  1 


Neural  Network  Response  vs.  Actual  Response  for  a0=2,b0=0.1,Rho=0.0001 
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Figure  4.8:  WNN  Design  6  -  Testing  Run  2 


51 


Neural  Network  Response  vs.  Actual  Response  for  a0=2,b0=0.1  ,Rho=0.0001 


time  (seconds) 


time  (seconds) 


Figure  4.9:  WNN  Design  6  -  Testing  Run  4 
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5,  MODELING  AND  CONTROL  OF  STRUCTURAL  SYSTEMS 

USING  WAVELETS 


In  this  chapter,  we  explore  a  couple  of  uses  for  time-frequency  wavelets.  These  include 
the  modeling  and  control  of  structural  systems  and  subsystems.  We  develop  a  versatile 
structure  model  based  on  time-frequency  wavelets  as  opposed  to  the  traditional  finite 
element/modal  models.  The  versatility  of  these  wavelet  models  provides  a  means  of 
developing  a  best  basis  for  a  given  structural  system.  These  wavelet  models  can  be 
easily  modified  to  account  for  structural  changes  due  to  the  operating  environment, 
coupled  dynamic  loading  and  material  degradation.  The  wavelet  models  are  developed 
for  the  purpose  of  being  utilized  during  the  controller  design  phase.  We  then  present 
some  theoretical  results  necessary  for  developing  a  sensitivity  minimization  controller 
based  on  a  wavelet  model.  The  theoretical  work  for  other  controller  objectives  can 
also  be  performed  as  extensions  of  this  work.  Initial  proof  of  concept  studies  will  be 
performed  on  several  flexible  structure  models  using  simulation. 

The  main  concept  presented  in  this  chapter  is  to  use  a  wavelet  basis  to  model  the 
structural  system  and  then  determine  the  optimal  controller  using  this  same  wavelet 
basis  in  an  effort  to  obtain  a  nesting  property  between  controllers  of  different  order. 
The  nesting  property  which  wavelets  possess  is  one  of  their  primary  strengths.  The 
control  problem  to  be  addressed  is  that  of  minimizing  the  effect  of  disturbances  on 
a  SISO  closed  loop  system,  i.e.,  minimizing  the  sensitivity  function,  in  terms  of  the 
norm.  This  control  problem  solution  will  be  discussed  herein  and  appeared  in 
[46].  The  time/space  localization  property  of  wavelets  provides  for  the  capability 
of  “localized”  control  action  over  a  selected  period  of  time.  This  type  of  control 
should  prove  invaluable  for  structural  systems  with  stringent  positional  or  pointing 
reciuirements. 

Wavelets  provide  a  rich  theory  with  a  library  of  functions  and  algorithms  available 
for  modeling  signals.  Wavelets  also  have  the  al^ility  to  attain  a  best  basis  for  a 
signal  or  class  of  signals.  This  means  that  a  best  basis  of  wavelet  functions  can  be 
easily  modified  to  account  for  uncertainties,  Failures,  and  changes  in  the  operating 
environment  of  a  structural  system  by  simply  changing  the  coefficients  on  the  wavelet 
basis  functions.  Since  wavelet  coefficients  tend  to  decay  ciuickly  and  usually  recpiire 
substantially  fewer  basis  functions  compared  with  a  Fourier  basis,  truncation  of  an 
infinite  dimensional  wavelet  model  to  a  finite  diinensional  model  should  introduce 
negligible  modeling  error.  In  addition,  if  the  controller  also  depended  on  the  same 
wavelet  basis  functions  as  the  structural  system,  then  the  controller  could  also  be 
easily  modified  to  account  for  uncertainties,  changes  in  operating  environment  of  the 
structure,  etc.,  making  for  a  more  robust  controller. 
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5.1 


Background  on  Time-Frequency  Wavelets 


There  are  many  excellent  articles  and  books  on  general  wavelet  theory,  for  example, 
see  [37],  [47],  [48],  [49],  [40],  [50]  and  [41]  and  the  references  therein.  Wavelets  are  a 
set  of  basis  functions,  tpm.ni^)  which  can  be  used  to  represent  all  admissible  functions 
f(x)  G  H,  Hilbert  space,  as  follows: 

fix)  =  (ST) 

m,n 


The  interesting  characteristic  of  wavelets  is  that  the  basis  functions,  0m, con¬ 
structed  from  dilations  and  translations  of  a  mother  wavelet,  0(x).  Time- frequency 
wavelets  are  closely  linked  to  the  Short  Time  Fourier  Transform  (STFT)  which  was  pi¬ 
oneered  by  Gabor  [51],  However,  ‘‘Gabor”  wavelets  suffer  from  serious  algorithmic  dif¬ 
ficulties  and  can  never  be  an  orthonormal  basis  for  C^{1R.).  A  significant  improvement 
over  the  “Gabor”  wavelets  has  been  the  development  of  orthogonal  time-frequency 
wavelets  by  Coifman  and  Meyer  [52],  and  Malvar  [53,  54].  These  researchers  use 
smooth  window  functions  to  split  the  signal  and  employ  extension  concepts  from  the 
discrete  cosine  transform  (DCT-IV)  to  “fold”  the  overlapping  parts  back  into  the  split 
signal.  Their  clever  construction  allows  one  to  construct  smooth  orthonormal  bases 
for  C^{1R)  which  depend  on  arbitrary  partitions  of  the  real  line  [55],  [48],  [49]. 

The  following  definitions  are  essential  to  understanding  the  construction  of  orthonor¬ 
mal  time-frequency  wavelets.  Some  of  these  definitions  can  be  found  in  Wickerhauser, 
[48]. 

Definition  2  (Function  Polarity)  A  function  f  has  positive  polarity  (+)  loith  re¬ 
spect  to  the  point  a  on  the  interval  [a  —  (.a  +  c]  f{t)  ~  /(2a  --  t)  for 

t  G  [a  —  e^a  +  e].  In  other  words,  f  is  an  even  function  about  the  point  a.  A  function  f 
has  negative  polarity  (-)  with  respect  to  the  point  a  on  the  interval  [a  —  e,  a  +  e] 

/(f)  =:  — /(2a  — f)  fort  G  [a  — e,a  +  e].  In  other  words,  f  is  an  odd  function  about  the 
point  a. 


Definition  3  (Cutoff  Functions)  A  cutojf  fuiiclion,  m{t)  6  is  a  complex, 

continuous,  d-times  differentiable  function  with  the  following  properties: 


m{t)\^  +  |m(— =  1  'it  ^  R 


m{t)  = 


0,  if  t  <  -1, 
1,  if  t  >  1. 
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Definition  4  (Bell  Function)  A  bell  function,  bj{t),  supported  on  the  interval  [oj  — 
ej,aj+i  +  Cj+i]  is  given  by: 

bj{t)  —  <  1,  i  €  [«i  +  ej,aj+i  —  Cj+i] 

.  "^i+i  ^  ^  “j+i  +  ej+i] 


Definition  5  (Folding  Operator)  The  folding  operator,  F ,  is  given  by: 

'  -  t),  if  Gj  -Cj  <t<  aj 

F{m)g{t)  =  <  m{^-^)g{t)  +  m{^)g{2aj  -  t),  if  aj  <t  <  aj  +  ej, 

g{t)^  otherwise 

Definition  6  (Unfolding  Operator)  The  adjoint  of  the  folding  operator^  F*  the 
unfolding  operator^  is  given  by: 

[  +  ^(^)5(2aj  -t),  if  aj  -tj<t<  aj 

F*{m)g{t)  =  I  m{^-^)g{t)  -  m{^)g{2aj  -  t),  if  aj  <  t  <  aj  +  e^-, 

[  g{t),  otherwise 

Unitary  Operator  Property:  F*{m)F{ni)g[t)  ~  F{m)F*{m)g{t)  Vt  ^  aj 
Action  Interval:  [aj  —  Cj ,  aj  +  ef\ 

Definition  7  (Translation  Operator)  The  translation  operator  is  given  by: 
SMt)=s(t-a) 

The  adjoint  of  the  translation  operator  is  given  by: 

Slg{t)  =  g{t  +  a) 

Definition  8  (Dilation  Operator)  The  dilation  operator  is  given  by:  Deg{t)  = 
The  adjoint  of  the  dilation  operator  is  given  by:  D*^g[t)  =  e^l'^g[et) 


The  translation  and  dilation  operators  can  now  be  used  to  make  the  action  region  of 
the  folding  and  unfolding  operators  the  interval  (a-e,a+e)  as  follows: 

F{mk,ak,ek)g{t)  =  Sa,De,F  Dl^Sl^g{t) 

\  hik{^)g{t)  -  mkC-^)g{2ak  -  t),  if  -  Ck  <  t  <  Ok 
=  j  ^rikC-^)g{t)  +  mk{^)g{2ak  -  t),  ifak<t<ak  +  ek,  (5.2) 

(  g{t),  otherwise. 
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(5.3) 


F"{mk,ak,ek)g{t)  =  Sa,De,F*  Dl^S*^g{t) 

(  i^k{^)g{t)  +  mkC-^)g{2ak  -  i),  if  ai,  -  ek  <  t  <  cik 
=  \  rnkC-^)g{t)  -  m{^)g[2ak  -  t),  if  cik  <  t  <  ctk  +  tk, 
[  g[i)^  otherwise. 


Notation:  Let  F{mk,ak,ek)  =  Fk  and  F*{mk,ak,ek)  =  F^ 

Definition  9  (Disjoint  Folding  Operators)  A  family  of  folding  {Fk}  and  unfold¬ 
ing  {Fj^}  operators  is  disjoint,  and  therefore  commutes,  if 
{ak  O-k,  Ok  “I"  ^A:)r*l(^A:4'l  5  ^A:+l  ~f” 


Now,  the  real  line,  IR,  is  partitioned  into  compatible,  disjoint  intervals,  M.  =  U^-  h 
and  PIa:  h  =  0  as  follows: 

Let  {ck  :  k  E  Z}  he  any  sequence  such  that 


•  For  k  <  j ,  Ck  <  Cj 

•  Ck  oo  and  C-k  — oo  as  A:  — >  oo 

For  each  k  ^  Z,  define: 


-[cfc  +  CA+l] 

(5.4) 

-[Cfc+I  -  Ck\ 

(5.5) 

[ak,ak+i] 

(5.6) 

Figure  5.1  shows  how  the  intervals  for  the  trigonometric  basis  are  partitioned  and 
demonstrates  that  any  given  interval  Ik  is  only  overlapped  by  the  intervals  Ik-i  and 
Ik+i-  Figure  5.1  also  demonstrates  that  the  midpoint  of  the  ^h  interval,  Ik,  is  Ck+i- 

This  partition  can  then  be  refined  by  splitting  the  intervals  at  their  midpoints.  For 
the  jth  refinement  level,  the  intervals  are  given  by: 


ak  ak+i 
21 


(5.7) 


The  development  of  local  trigonometric  transforms  which  serve  as  a  complete  or¬ 
thonormal  basis  for  C^{IR)  hinges  on  the  properties  of  the  bell  functions,  the  def¬ 
initions  of  the  folding  and  unfolding  operators,  and  the  polarity  of  trigonometric 
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Adjacenl  Intervals  With  Disjoint  Action  Regions 


<- 


- 


l(k) 


IIk+1] 

- 


Figure  5.1:  Partition  of  intervals,  {/a:},  constructed  using  disjoint  action  regions  at 
the  interval  endpoints 

transform  used.  Two  examples  of  complete,  orthonormal  local  trigonometric  bases 
are  the  Local  Cosine  Transform  of  Type  IV  (LCT-IV)  and  the  Local  Sine  Trans¬ 
form  of  Type  IV  (LST-IV).  When  the  bells  bj{t)  are  symmetric  with  respect  to  their 
centers,  i.e.,  ej  =  Vj,  and  the  same  cutoff  function  is  used  throughout,  i.e., 

mj{t)  =  Vj,  the  interval  I  =  (a^^,  ak^i)  is  compatible  with  its  adjacent  trans¬ 
late  by  |/|  =  —  a^.  For  this  special  case,  the  folding  and  unfolding  operators 

simplify  as  follows: 


Definition  10  (Folding  Operator  for  Compatible  Intervals)  The  folding  oper¬ 
ator  for  compatible  intervals^  FC,  is  given  by: 

aj+i -t),  if  a,  <  t  <  aj  +  e, 

+  ttj+i  -  t),  if  aj+i  -  e<t  <  aj+i 
g{t),  otherwise 


FC{m)g{t)  =  \ 


Definition  11  (Unfolding  Operator  for  Compatible  Intervals)  The  unfolding 
operator  for  compatible  intervals,  FC",  is  given  by: 

m{^-^)g{t)  -  m{^)g{aj  +  Uj+i  -  t),  if  aj  <  t  <  aj  +  e, 

mC^^^-)git)  +  rhC-^-^)9ich  +  «j+i  “  if  «j+i  -  e  <t  <  a,+i 

g{t),  otherwise 


FC*{m)g{t)  =  { 


With  the  use  of  these  compatible  interval  folding  and  unfolding  operators,  there  is  no 
longer  any  need  to  step  outside  of  the  interval  when  performing  folding  and  unfolding. 

The  general  form  of  a  typical  family  of  time-frequency  wavelets  is  given  by: 


57 


fY'  /tt  (n  +  i)  (<  -  ftfc)  , 

=  Fi^^k+lXM  COSn,k{t)  =  ^ I  - [7^1 -  I 


where:  4\^k{t)  G  L^{M) 

xM  =  xC-^)  = 


1,  if  (Ik  <  t  <  Uk+U 
0,  otherwise. 

^(n+|)(i-at) 


cos„..(0  = 

Wk{t)  =  mk  (^)  ^k+i  (^^),  Supp  iUk{t)  =  [uk  -  ek,ak+i  +  ek+i] 
h  =  [o-k^<^k+i]  is  the  kth.  interval 
\h\  =  Ofc+i  —  «/; 


The  time-frequency  wavelets  allow  for  windows  of  varying  length  in  which  the  time 
domain  is  divided  into  unequal  intervals.  Figure  5.2  shows  a  typical  time-frequency 
wavelet. 


Local  Cosino  Wavelal 


Figure  5.2:  Typical  Time- Frequency  Wavelet 
Properties  of  f>n,k{t)-. . 

(i)  Orthonormality:  {ifn,k{t)-,‘4’m,j{t))  =  8„,,nhk.j- 

(ii)  Compact  Support:  Supp  ipn.kit)  = 

[u/c  ^k^^k+l  T  6fc-|-l]  "  ''  V’n,A:(t)  G  L  (IR). 

(iii)  Minimal  Overlap:  For  a  given  interval,  4,  i.'n.k{t)  is  only  overlapped  by  ipn.k-iit) 
and  fin,k+\{t)  due  to  being  constructed  from  disjoint  folding  and  unfolding  op¬ 
erators. 
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Remark 

The  construction  of  disjoint  intervals  for  ^(0,oo)  requires  that  c_^  — ^  0  as  k 
oo.  Such  a  partition  yields  a  continuous,  orthonormal  local  trigonometric  basis  for 
L^(0,oo)  [55].  Now.  for  appropriately  selected  cutoff  functions,  an  orthonormal 
wavelet  basis  can  be  determined  by  taking  the  Laplace  transform  of  each  element  in 
this  basis  for  T^(0,oo). 


5,2  Structural  System  Modeling  and  Identification 

The  control  of  advanced  structural  systems  is  a  demanding  task  due  mainly  to  the 
intrinsic  dynamic  characteristics  of  these  structures.  Typical  modeling  approaches  in¬ 
volve  using  structural  analysis  programs  like  NASTRAN  for  developing  finite  element 
models  [56,  57].  Typically,  finite  element  models  can  be  characterized  by  the  presence 
of  densely  packed  modal  frequencies  with  low  damping  [58].  In  addition,  finite  element 
models  do  not  account  for  the  many  uncertainties  present  in  the  structure,  including 
structural  uncertainties,  failure  of  system  components,  and  changes  in  the  operating 
environment.  For  a  1%  error  in  the  dominant  modal  frequency,  the  modeling  error 
could  be  as  large  as  the  contribution  of  the  mode  itself  [59].  Truncation  and  round 
off  inherent  in  using  a  finite-dimensional  model  for  an  infinite  dimensional  system 
invariable  introduce  additional  errors  in  the  modal  frequencies  and  the  mode  shapes. 
Lastly,  different  modes  may  be  excited  during  different  operating/failure  conditions 
requiring  a  versatile  model  of  the  structure. 

A  general  SISO  structural  plant  model  can  be  represented  as: 


where  gk  is  a  gain  term.  Assume  that  any  rigid  body  modes  have  been  stabilized  by 
an  inner  loop.  Then  P(s)  G  The  structural  plant,  P(s),  can  now  be  modeled  by 
using  a  time-frequency  wavelet  representation  of  the  form: 


P{s)  =  Y,Otk'^n,,k{s)  (5.10) 

where  ^rik,k{^)  is  an  wavelet  basis  funct  ion.  Due  to  the  linearity  of  the  Fourier 
Transform,  and  the  fact  that  the  Fourier  'rraiisforin  is  an  isomorphism  between 
Z/^[0,oc)  and  7/^,  the  equivalent  time  domain  [('presentation  of  structural  system 
is  given  by: 


59 


(5.11) 


OO 

P(^)  = 

k=0 

where  i>nk,k{t)  €  L'^{0,oo)f]L^{0,oo). 

The  control  designer  would  like  to  be  able  to  obtain  a  best  basis  for  a  given  structural 
system,  or  possibly  a  best  basis  (BB)  for  a  class  of  structural  systems  /  subsystems. 
A  best  basis  selection  algorithm  [60],  based  on  a  signal  processing  viewpoint,  has 
already  been  developed.  The  BB  selection  procedure  is  an  adaptive  algorithm  which 
determines  the  splitting  locations  for  the  time-frequency  wavelets  based  on  the  signal 
to  be  decomposed.  The  decision  is  based  on  yielding  a  representation  which  has  the 
most  decorrelation.  Entropy  is  one  method  out  of  many  possible  methods,  [48],  used 
to  measure  the  decorrelation.  There  are  also  many  other  adaptive  techniques  which 
can  be  utilized  in  the  selection  of  the  best  basis  for  a  given  signal.  It  should  be  noted 
that  for  the  controller  design  presented  herein,  the  BB  selection  procedure  would 
have  to  be  modified  to  ensure  that  the  basis  elements  selected  remain  orthogonal  to 
each  other.  Another  possible  approach  is  to  determine  the  best  level  basis  out  of 
the  total  number  of  depth  levels  considered  for  the  given  structural  system.  This 
best  level  basis  could  serve  as  the  structural  basis  and  would  maintain  the  recpiired 
orthogonality  between  the  basis  elements.  Figure  5.3  provides  an  example  of  possible 
splitting  locations  for  a  time-frequency  wavelet  representation  of  a  particular  signal 
and  then  gives  a  potential  orthogonal  basis  and  a  potential  level  basis,  which  by 
definition  is  guaranteed  to  be  orthogonal. 

Figure  5.4  outlines  the  procedure  for  the  decomposition  of  a  particular  discrete  signal, 
x{k),  into  cosine  packets  (Cosine  Packet  Analysis)  and  the  procedure  for  the  recon¬ 
struction  of  the  discrete  signal  from  the  cosine  packets  (Cosine  Packet  Synthesis). 
Figure  5.4  demonstrates  these  procedures  for  a  given  depth  level,  d,  in  the  packet  ta¬ 
ble,  corresponding  to  a  level  basis  analysis/synthesis  procedure.  Note  that  DCT-IV 
denotes  the  discrete  cosine  transform  of  type  IV.  It  should  be  noted  that  sine  packets, 
using  the  discrete  sine  transform  of  type  IV  (DST-IV),  or  a  combination  of  both  sine 
and  cosine  packets  can  also  be  used. 

Figures  5.5  -  5.7  presents  an  example  of  a  time-frequency  wavelet  decomposition, 
in  particular,  the  cosine  packet  decomposition,  basis  selection,  and  cosine  packet 
reconstruction  of  the  impulse  response  for  a  structural  system  containing  two  modes. 
The  basis  to  be  determined  was  the  best  basis  for  the  system  with  entropy  as  the 
information  cost  function  to  be  minimized  by  the  basis  selection  procedure. 

Figures  5.8  -  5.10  presents  the  same  example  of  the  cosine  packet  decomposition,  etc., 
as  figures  5.5  -  5.7,  but  the  basis  to  be  determined  is  the  best  level  basis  for  the  system 
with  entropy  as  the  information  cost  function. 
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5.3 


Multiresolutional  Sensitivity  Minimization 


The  feedback  control  system  is  shown  in  Figure  5.11,  where  P  is  the  plant  to  be  con¬ 
trolled  and  W  is  the  weight  modeling  the  disturbances.  The  Wiener-Hopf  approach 
to  the  sensitivity  minimization  problem  involves  finding  an  internally  stabilizing  con¬ 
troller,  C,  such  that  the  following  performance  is  achieved: 


inf.  py5||, 

G  staotlizing  r 


The  sensitivity  function  for  the  feedback  control  system  is  given  by: 
S  :=  {l  +  PC)-^ 


(5.12) 


(5.13) 


The  sensitivity  minimization  problem  has  been  studied  by  numerous  researchers  in 
terms  of  minimizing  the  norm  and  the  H'^  norm,  see  [61],  [62]  and  [6.3]  and  the 
references  contained  therein. 

Assumptions 


1.  We  RH\  i.e.,  W  is  real,  rational,  stable  and  strictly  proper. 

2.  P  e  H^,  i.e.,  P  is  stable. 


The  first  step  in  our  MC  development  is  to  perform  a  Youla  parameterization  [64]  for 


the  controller  C  according  to: 

C  =  Qil-PQ)-\  QeH^,  {l-PQ)^0 

(5.14) 

The  performance  measure,  equation  (5.12),  becomes 

^  11^(1 

(5.15) 

=  inf  \\W-PQ)\\2 

(.5.16) 

where:  Q  =:WQ 

Now  represent  the  performance  criteria,  equation  (5.16), 

as  an  binner  product: 

J  =\\W  -  PQ\\l  =  {W  -  PQ,W  -  PQ) 

=  {W,  W)  -  2Re{{W,  PQ)]  4-  (PQ,  PQ) 

(5.17) 
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Let  the  Fourier  Transform  of  a  signal  be  represented  by  X{juj)  and  recall  the 
following  theorem: 

Theorem  8  (Plancherel)  If  fit)  G  L^iM),  then  Fijtn)  G  and  ||/||  =  \\F\\. 

In  addition,  {f,g)  =  {F,G)  for  any  two  functions  f,g  G  L^{IR). 

In  light  of  this  theorem,  the  performance  criteria,  J,  can  be  represented  equivalently 
as: 


J  =  {w,w)  -  2Re{fw,pq)]  +  {p(l,pq)  (5.18) 

The  next  step  in  our  MC  development  is  to  represent  the  plant  P  using  an  orthonor¬ 
mal,  preferably  best,  wavelet  basis  according  to  equation  (5.10)  given  previously. 
Now,  the  disturbance  weight  W,  and  the  controller  function  Q  will  also  be  represented 
by  the  same  orthonormal,  wavelet  basis  as  follows: 


W{s)  =  Eidk^n.As)  (5.19) 

^'=:0 

0(5)  =  f:7/.-'J..,/.-(5)  (5.20) 

k=0 

The  equivalent  time  domain  representation  of  these  signals  is  given  by: 

<-‘^it)  =  ^  (5.21) 

k=Q 

q{t)  =  (5-22) 

k=0 

Define: 

r  [7o7i72  •  • •  (5.23) 


The  performance  measure  is  to  be  minimized  by  choosing  the  optimal  gains,  P,  for  the 
controller  function  Q.  However,  due  to  the  local  nature  of  the  wavelet  basis  functions 
used  to  represent  the  components  of  the  system,  the  performance  measure,  eciuation 
(5.18),  can  be  written  as: 


./(r)  =  ^j,(7fc) 

4 


(5.24) 
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Making  use  of  equation  (5.24)  and  Figure  5.1  of  Section  5.1,  the  optimal  controller 
gains,  are  determined.  The  complete  derivation  is  presented  in  [46].  The  optimal 
controller  gain  for  a  given  interval,  k,  can  be  represented  functionally  as: 

Ik  =  li'^k-l,0(k,ak+l,  Plc-1,  l3k,l^k+l,4’k-l,i''k,4’k+U'fk-\,"fk+l)  (5.25) 

5.4  Multiresolutional  Controller  Development  and  Implementation 

The  theoretical  development  of  a  multiresolutional  controller  for  sensitivity  minimiza¬ 
tion  was  presented  in  Section  5.3.  The  section  discusses  some  of  the  implementation 
issues  prevalent  to  the  design  of  a  multiresolutional  controller  (MC)  for  sensitivity 
minimization.  First,  the  functional  representation  of  the  controller  gains,  jk,  found 
in  equation  (5.25),  actually  represents  inner  product  terms  of  the  time- frequency 
wavelets,  ^/’fc(t),  which  can  be  evaluated  numerically  using  equation  (5.8)  and  the 
unitary  operator  property  of  the  folding  operators,  see  Definition  5,  as  follows: 

=  {Fk+iFkf,Xh{t)(^oSn,k{t))  (■5-26) 

Secondly,  the  optimal  controller  gains  are  a  function  of  the  best  basis  chosen  for  the 
plant.  Therefore,  the  controller  design  is  extremely  problem  and  plant  specific.  This 
is  by  no  means  a  drawback,  and  actually  is  exactly  what  a  control  designer  would 
want.  However,  a  suboptimal  controller  should  be  directly  related  to  the  optimal 
controller  when  represented  using  a  wavelet  basis.  Determining  the  optimal  algorithm 
to  use  for  the  selection  of  the  best  wavelet  basis  functions  also  needs  to  be  addressed. 
However,  the  best  basis  selection  algorithm  [60]  appears  to  be  a  good  starting  place, 
keeping  in  mind  that  the  resultant  best  basis  must  be  orthogonal  so  that  the  controller 
methodology  developed  is  applicable.  Next  is  the  need  to  truncate  the  infinite  wavelet 
basis  representation  for.  the  plant.  Again,  this  is  not  a  drawback,  but  a  feature  of 
wavelets  that  is  usually  exploited,  especially  in  \'id('o  compression  problems.  Wavelets 
are  such  a  powerful  tool  since  they  require  loss  basis  functions  than  the  standard 
Fourier  methods  due  to  their  localization  in  both  time  and  frequency.  Another  issue 
is  selection  of  the  best  optimization  technique  to  use  when  solving  equation  (5.25)  for 
the  optimal  7jt’s.  Due  to  the  coupling  between  the  and  terms,  a  recursive 

solution  is  possible  which  would  relate  the  first  gain  to  the  last  7  gain,  and  then 
all  the  intermediate  gains  would  be  found  by  back  substitution.  Other  optimization 
procedures  could  also  be  employed,  possibly  a  dynamic  programming  method  starting 
at  the  last  time  interval  or  a  conjugate  gradient  |n‘ocedure.  Lastly,  suboptimal  7a:’s 
could  always  be  calculated  by  simply  calculating  bounds  on  the  cross  coupling  terms, 
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and  using  these  bounds  in  the  calculation  for  the  7t’s.  Then  the  optimization  problem 
is  directly  solved  by  determining  the  7^-  gains,  for  all  k  =  0,1,  N,  using  equation 
(5.25). 
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Potential  splitting  locations 
for  a  time-frequency 
wavelet  representation 
over  an  interval  I 

Depth  Level _ 

d  =  0  _ _ 

d  =  l 


d  =  3  _  - 

- I 


Orthogonal  Basis 


d  =  0 
d  =  l 
d  =  2 
d  =  3 


Level  Basis 


d  =  0 
d  =  l 
d  =  2 
d  =  3 


Figure  5.3:  Example  of  possible  splitting  locations  for  a  particular  signal  and  a  po¬ 
tential  orthogonal  basis  and  a  potential  level  basis 
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COSINE  PACKET  ANALYSIS:  For  a  given  depth  level,  d 


x(k) 


For  b  =  0, 2 
Partition  into 
Compatible 
Packet  Intervals 


j  I[b-1] 

I _ IW_ 

i  I[b+1] 


►j  DCT-IV 


CosPkt|j(n) 


COSINE  PACKET  SYNTHESIS:  For  a  given  depth,  d 
For  b  =  0, ...,  2 

;  I[b-1] 

p^T-IV  i  FC*(mj)  I[b]  ^*^2: 

‘  I[b+1]  < 


Figure  5.4;  Procedures  for  performing  cosine  packet  analysis  and  synthesis  on  a  dis¬ 
crete  signal,  x(k),  for  a  particular  depth  level,  d. 


Impulse  Response  for  Two  Modes  with  zeta  =  0.01  &  omega  =  6.283  &  12.57 


Frequency  Response  for  Two  Modes  with  zeta  =  0.01  &  omega  =  6.283  &  12.57 


Figure  5.5;  Impulse  and  frequency  response  for  a  fle.xible  structure  with  two  modes 
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Time(Frequency) 


Cosine  Packet  Tree  for  the  Best  Basis  Time-Frequency  Plane  for  the  Best  Basis 


Figure  5.6:  Cosine  packet  decomposition  and  best  basis  selection  for  a  flexible  struc¬ 
ture  with  two  modes  at  =  27r  and  (jJ2  =  47r  with  the  identical  damping  factor  of 

C  =  o.oi. 


67 


Impulse  Response  for  Two  Modes  with  zeta  =  0.01  &  omega  =  6.283  &  1 2.57 


.jq-14  Reconstruction  Error  of  the  Impulse  Response 


Figure  5.7:  Reconstruction  of  the  impulse  response  for  a  flexible  structure  with  two 
modes 


Impulse  Response  for  Two  Modes  with  zeta  =  0.01  &  omega  =  6.283  &  12.57 


30  40  50 

Time  in  seconds 


Frequency  Response  for  Two  Modes  with  zeta  =  0.01  &  omega  =  6.283  &  12.57 


Figure  5.8:  Impulse  and  frequency  response  for  a  flexible  structure  with  two  modes 
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0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

Time(Frequency) 


Cosine  Packet  Tree  for  the  Best  Level  Basis  Time-Frequency  Plane  for  the  Best  Level  Basis 


Figure  5.9:  Cosine  packet  decomposition  and  best  level  basis  selection  for  a  flexible 
structure  with  two  modes  at  lji  =  27r  and  uj2  =  47r  with  the  identical  damping  factor 
of  C  =  0.01. 
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Impulse  Response  for  Two  Modes  with  zeta  =  0.01  &  omega  =  6.283  &  12.57 

20 
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0  10  20  30  40  50  60  70  80 

Time  in  Seconds 

Reconstructed  Signal  from  the  Best  Level  Basis  Cosine  Packets 
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’^^0  10  20  30  40  50  60  70  80 

Time  in  Seconds 

j.  .jQ-’4  Reconstruction  Error  of  the  Impulse  Response 


Figure  5.10:  Reconstruction  of  the  impulse  response  for  a  flexible  structure  with  two 
modes 


Figure  5.11:  Feedback  Control  System 
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6.  CONCLUSION 


The  main  thrust  of  this  project  has  been  to  consider  the  control  of  interconnected  rigid 
and  flexible  structures.  Due  to  the  fact  that  we  distributed  our  attention  to  different 
modeling,  analysis  and  design  techniques,  we  have  obtained  many  new  results  in 
different  areas. 

Conclusions  obtained  from  individual  aspects  of  our  report  are  listed  below,  as  they 
correspond  to  the  chapters  in  this  report: 

In  Chapter  2  several  results  for  sliding  mode  control  of  sampled-data  systems  were 
outlined.  Due  to  limited  sampling  frequency  and  the  ‘fliold”  processes  in  the  control 
channels,  perfect  sliding  mode  never  occurs  in  a  sampled-data  system.  The  control 
law  being  piecewise  constant,  the  states  can  only  be  kept  in  a  boundary  layer  of 
the  sliding  surface  instead  of  being  constantly  and  exactly  on  that  surface.  In  this 
chapter,  a  new  method  of  implementing  discrete  time  sliding  mode  control  for  sample- 
data  systems  was  developed  to  reduce  the  thickness  of  the  boundary  layer,  to  alleviate 
the  effect  of  uncertainties,  and  to  remove  the  undesired  chattering  in  sliding  mode. 
Three  classes  of  sampled-data  systems  with  uncertainties  were  investigated;  namely, 
linear,  nonlinear,  and  stochastic  systems.  By  proper  consideration  of  the  sampling 
phenomenon  in  the  control  design,  one  can  maintain  the  states  in  the  vicinity  of  the 
sliding  surface  with  accuracy  up  to  at  least  0{T^).  and  0{T^)  respectively. 

We  have  demonstrated  these  results  on  an  actual  hardware  structure. 

In  Chapter  3,  we  have  introduced  a  new  approach  for  vibration  damping  in  flex¬ 
ible  structures  based  on  the  sliding-mode  control  approach.  Previous  use  of  this 
approach  has  been  for  finite  dimensional  (approximate)  models  of  flexible  structures. 
Here,  we  retain  the  infinite-dimensional  model  for  the  structures  and  investigate  exact 
solutions  using  sliding-mode  control.  In  the  particular  strategy  introduced  here,  the 
design  is  accomplished  initially  for  a  class  of  differential-difference  systems.  It  is  then 
demonstrated  that  the  model  of  flexible  structures  with  second  order  spatial  partial 
derivatives,  can  be  transformed  into  the  above  differential-difference  form.  A  number 
of  different  cases  have  been  analyzed.  We  then  consider  a  class  of  flexible  structures 
with  fourth  order  spatial  partial  derivatives.  An  integral  transform  is  introduced  that 
changes  the  model  into  second  order  form.  Thus,  a  two-stage  design  process  can  be 
utilized  to  eventually  generate  the  sliding-mode  controllers.  Sliding  mode  algorithms 
considered  make  the  closed  loop  system  highly  insensitive  to  external  disturbances 
and  parameter  variations.  Furthermore,  from  the  applications  viewpoint,  the  models 
considered  are  particularly  appropriate  for  the  ut  ilization  of  distributed  actuation  and 
may  be  used  in  ^^smart  structures”  with  piezo(dectric  materials.  To  avoid  difficulties 
in  solving  for  the  transformation  kernel  and  performing  the  integral  transformation 
analytically,  a  neural  network  may  be  trained  to  perform  these  tasks.  By  using  the 
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neural  network,  exact  knowledge  of  the  kernel  of  the  transformation  is  not  required. 
The  neural  network  then  becomes  the  first  stage  of  the  vibration  suppression  con¬ 
troller.  The  second  stage  of  the  composite  controller  is  the  sliding  mode  controller 
which  has  been  designed  for  the  differential-difference  system. 

In  Chapter  4,  we  introduce  a  new  general  design  procedure  for  radial  basis  function 
neural  networks  and  wavelet  neural  networks.  These  type  of  neural  networks  are 
especially  well-suited  for  solving  general  pragmatic  learning  problems  in  which  the 
training  data  consists  of  irregularly  positioned  data  points.  A  multiresolution  learn¬ 
ing  network  or  wavelet  network,  has  the  added  capacity  of  being  able  to  generalize 
the  nonlinear  relationship  between  the  given  inputs  and  desired  outputs,  along  with 
capturing  the  desired  fine  details  of  the  nonlinear  mapping.  VVe  present  simulation 
results  for  these  two  types  of  networks  based  on  predicting  actuator  nonlinearities  oc¬ 
curring  in  a  flexible  structure.  The  neural  networks  were  trained  using  a  set  of  training 
data  and  were  then  tested  using  additional  data  not  contained  in  the  training  set  and 
excellent  performance  was  achieved  and  demonstrated. 

In  Chapter  5,  we  explored  a  couple  of  uses  for  time-frequency  wavelets.  First, 
an  overview  of  time-frequency  wavelets  was  presented.  We  then  presented  the  con¬ 
cept  of  modeling  a  structural  system  using  a  best  basis  made  up  of  time-frequency 
wavelets.  We  developed  a  versatile  structure  model  based  on  time-frequency  wavelets 
as  opposed  to  the  traditional  finite  element/modal  models.  The  decomposition  and 
reconstruction  schemes  for  time-frequency  wavelets  were  presented.  Then  simulation 
results  were  presented  for  modeling  the  impulse  response  of  a  structural  system  using 
cosine  packets.  These  simulation  results  demonstrated  the  cosine  packet  decompo¬ 
sition,  the  basis  selection  procedure,  and  the  reconstruction  process.  The  concept 
of  a  multiresolutional  controller  was  formulated  and  discussed.  A  multiresolutional 
controller  which  minimizes  the  sensitivity  function  for  a  SISO  flexible  structure  in  the 
norm  was  derived.  This  controller  utilizes  the  structural  models  which  were  devel¬ 
oped  using  time-frequency  wavelets.  Implementation  issues  for  the  multiresolutional 
controller  were  also  discussed. 
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